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^ ' Abstract. We show that almost all the linear differential operators factors 

1/^ ' obtained in the analysis of the n-particle contributions x'"''s of the susceptibility 

ff^ ^ of the Ising model for n < 6, are linear differential operators "associated 

l/~\ , with elliptic curves". Beyond the simplest differential operators factors which 

f*"^ ■ are homomorphic to symmetric powers of the second order operator associated 

• ' with the complete elliptic integral E, the second and third order differential 

^^ ' operators Z2, F2, -F3, -£3 can actually be interpreted as modular forms of 

^ ^ ' the elliptic curve of the Ising model. A last order-four globally nilpotent linear 

CD , differential operator is not reducible to this elliptic curve, modular forms scheme. 

This operator is shown to actually correspond to a natural generalization of 
this elliptic curve, modular forms scheme, with the emergence of a Calabi- 
Yau equation, corresponding to a selected 4F3 hypergeometric function. This 
k^ " hypergeometric function can also be seen as a Hadamard product of the complete 

KN . elliptic integral K, with a remarkably simple algebraic pull-back (square root 

extension) , the corresponding Calabi-Yau fourth-order differential operator having 
a symplectic differential Galois group 5P(4, C). The mirror maps and higher 
order Schwarzian DDEs, associated with this Calabi-Yau ODE, present all the nice 
physical and mathematical ingredients we had with elliptic curves and modular 
forms, in particular an exact (isogenics) representation of the generators of the 
renormalization group, extending the modular group SL{2, Z,) to a GL{2, Z) 
symmetry group. 
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1. Introduction 

In a previous paper [I] some massive computer calculations have been performed on 
the susceptibility of the square Ising model and on the n-particle (n-fold integrals) 
contributions x^") of the susceptibility f2l|3l|4]. In three more recent papers [5l[6l[7j the 
linear differential operators for x and x^^' were carefully analysed. In particular, it 
was found that the minimal order linear differential operator for x'-^'' can be reduced to 
a minimal order linear differential operator L29 of order 29 for the linear combination 

This specific linear combination series being annihilated by an ODE of lower order, 
one has, thus, the occurrence of a direct sum structure. It was found [F, "B^ that this 
linear differential operator £29 can be factorised as a product of an order-five, an 
order-twelve, an order-one, and an order eleven linear differential operators 

-^29 — L^ ■ L12 ■ Li ■ Lii, (2) 

where the order-one linear differential operator Li has a rational solution and where 
the order-eleven linear differential operator has a direct-sum decomposition 

in = (Z2 • TVi) (SV2® (F3 • F2 • LI), (3) 

where Lf and A^i are order-one (globally nilpotenljjl) operators, Z2 is the second 
order operator occurring in the factorization of the linear differential operator [5] 
associated with x ^^d seen to corresponqjj to a modular form of weight one [5], V2 
is a second order operator equivalenl|§| to the second order operator associated with 
X*-^-* (or equivalently to the complete elliptic integral E), F2 and F3 are remarkable 
second order and third order globally nilpotent linear differential operators [5l[9]. 

The order-five linear differential operator L5 was shown to be equivalent to 
the symmetric fourth power of (the second order operator g|) Le corresponding to 
the complete elliptic integral E. These operators were actually obtained in exact 
arithmetics [S1[B]. The order-twelve linear differential operator L12 has been shown 
to be irreducible and has been proved to not be a symmetric product of differential 
operators of smaller orders (see [B] for details) . 

^ That is of the form ([TJ for order-one operators (see below). L^ has the simple rational solution 

«;2/(l-4«))2. 

f As well as the second order operator occurring in Apery's proof of the irrationality of C{3). 

§ In the sense of the equivalence of linear differential operators 1101 1111 (corresponding to the 

"Homomorphisms" command in Maple). We refer to this (classical) notion of equivalence of linear 

differential operators everywhere in the paper. In the literature the wording "operators of the same 

type" is also used |12l . 

I This second order operator does play a central role in the Gauss-Manin, or Picard-Fuchs, origin of 

the (sigma form of the) Painleve equations occurring for the two-point correlation functions of the 

Ising model [TallTil . 
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Similar calculations were actually achieved |7| on x*^^^ . For the linear combination 

$(6) ^ ^(6) _2^(4) +1^(2)^ (4) 

one obtains a linear differential operator of (minimal) order 46 which has the following 
factorization 

LiQ = Lq ■ L23 ■ il7, (5) 

with 

Lij = L^QLs® {L4,-L3-L2), (6) 

where D,j: denotesQ d/dx, and where the Li,n's {n — 1,2,3) are order-one linear 
differential operators (see Appendix A in [7]) and L2, £3 and Lg are [7] respectively 
equivalent (homomorphic) to i^;, to the symmetric square of Le and to the sym,m,etric 
fifth power of Le- 

The factorization (if any) of the order twenty-three linear differential operator 
L23 is beyond our current computational resources (see [7] for details). 

• Understanding the elementary factors: the modular form challenges 

Among these various globally nilpotent factors ^, of large order operators, 
one discovers order-one linear differential operators which, because they are globally 
nilpotent, are all of the form 

^ l_ dHRjx)) 

" n' dx ' ^' 

thus having A^-th root of rational functions solutions. One also discovers operators 
of various orders which are equivalent to symmetric powers of Le, the second order 
operator corresponding to complete elliptic integral of the first (or second kind), E 
or K, like V2 in ^, or L5 in ©, or L3 in ([B]), or ig in dS]), a remarkable second 
order operator Z2 having a modular form interpretation [9 , and a miscellaneous set 
of operators F2, F3, L3, L4, L12 and ^23- 

These last linear differential operators are not equivalent to symmetric powers of 
Le, and are still waiting for a modular form interpretation, if we think that all the 
globally nilpotent factors of these operators of the "Ising class" jTB] should have an 
interpretation in terms of elliptic curves in a modern sensqjj. These linear differential 
operators, emerging in the analysis of x , namely the second order operator F2 and 
the third order operator F^, are waiting for such a modular form interpretation, as 
well as L^ and L4 emerging in the analysis of x ■ The two remaining operators L12 
and L23 are too involved, and of too large order, for seeking for a possible modular 
form interpretation (up to equivalence). 

The purpose of this paper is to provide a mathematical interpretation of the 
F2, F3, L3, L4 elementary "bricks" of the n-fold integrals x*'^^ and x'"^\ in order to 
mathematically understand the Ising model. In this paper we will, as much as possible, 
use the same notations as in our previous papers [H [SJ [71 |F1 [5] . The linear differential 
operators, or equations, are in the x = w variable for the high temperature differential 

X In this paper we will use the notations Dt (resp. Dz) for d/dt (resp. d/dz). 

t Before Wiles' recent results, only elliptic curves with the property known as " complex 

multiplication" had been shown to be parametrised by modular functions (by Shimura in 1961). 
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operators and in the x = w^ variable for the fow temperature differential operators, 
where w = (1 + s'^)/s/2 and where s = sinh(2iir) with the standard notations for 
the Ising model. 

2. Modular form recalls 

2.1. Modular form, recalls : Z2 and Apery modular linear differential operator 

Let us introduce as in j9] the order-two Heunian operator which has the following 
solution Heun{8/9, 2/3, 1, 1, 1, 1; t): 

^9 /I 1 9\„ „ 3i-2 



.t t-l 9i-8/ {9t-8)(t-l)t' 

A simple change of variable (see equation (43) in |9J): 

—8 a; , , 

^ " {l~Ax){l-x)' ^' 

transforms H into the order-two linear differential operatoiyj: 

n, - Dl (9) 

1 + 10 x - 19 a;2 - 92 a;3 -h 12 a;^ -F 224 x^ - 64 x'^ 



6 



{1 + 3x + Ax^) {1 - 2x) {1 + 2x) {1 ^ 4:x) {1 ~ x) ■ X 
{l + 7x + Ax^){l-2x)^ 



■D^ 



(l-|-3x-F4a:2) {l-Axf{\-xf-x' 

We found ^9, that the second order linear differential operator Z2, occurring as 
a factor of the linear differential operator annihilating x'-'^' , is homomorphic |10[ ITT| 
to the operator ^. Recall (16) that the fundamental weight-1 modular forrn^ h^ for 
the modular group ro(Af) for A^ = 6, can be expressed as a simple Heun function, 
Heun{9/8, 3/4, 1, 1, 1, 1, — 1/8), or as a hypergeometric function: 

^J^ (10) 

((f + 6)-'' (t3 + 18 f2 + 84t + 24)3)1/12 V 7 

/r 1 5 , , , (t + 9f(t + 8)^t 

X 2F1 [-,-]; [1]; 1728 ^ ' ^ ^ ' 



'12' 12^'^ " (t -h 6)3 (i3-hl8i2-^84i-F 24)3 

which is solution of the order-two linear differential operator: 

, / 1 1 1 \ t+6 , , 

^* +(^+i + rr9)-^' +(rr8)(^T9)T7- ^^^^ 

simply related to H. Therefore, after some changes of variables, one can see the 
(selected) solution of Zi as a hypergeometric function (up to a Hauptmodul pull-back) 
corresponding to a weight-1 modular forn^ (namely /ig in [19 ). 

f Do note two minus sign misprints in the numerator of the Dx coefficient equation (44) of [9]: 

— lOx -(- IQa;-^ must be replaced by +10a; — 19a;^, see Jgjl. 

X The modular form he is also combinatorially significant: the perimeter generating function of the 

three-dimensional staircase polygons ^17^ can be expressed [9] in terms 1161 of /ig. The modular form 

ftg also occurs in Apery's study of C{3)- 

§ The simplest weight-1 modular form is 2-Pi{[l/12, 5/12], [1], J) = 12i/*l »?(t)2 J-1/12 where J is 

the Hauptmodul, r] is the Dcdckind eta function, and r is the ratio of periods (see (4.6) in [18]). 
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To sum-up 'Hx^ given by (jH]), has the foUowing solution: 

S= [^■M,y\ 2Fi([-L,l];[l]; X,), where: 

1728 x- (1 +3a; +4a;2)2(H-2 2;)6' ^ ' 

_ , ,,, X ■ (1 + 3 a; + 4 x^f (1 + 2 x)^ (1 - 4 xf (1 - xf 
^^ ^ ^^^^ (1 +7x + 4x2)3. P3 ' 

P = 1 + 237 a; + 1455 x^ + 4183 x^ + 5820 x^ + 3792 x^ + 64 x^ 

The solution of the operator Z2, in terms of hypergeometric functions, can now 
be understood from this hypergeometric function [up to a m,odular invariant pull-hack) 
structure. 

2.2. Recall on modular forms: Hauptmoduls and fundamental modular curve 

In several papers j20| I21j we underlined, for Yang-Baxter integrable models with 
a canonical genus-one parametrization [221 123] (elliptic functions of modulus k), 
that the exact generators of the renormalization group must necessarily identify with 
various isogenics 24' which amounts to multiplying, or dividing, r the ratio of the 
two periods of the elliptic curves, by an integer |25j . The simplest example is the 
Landen transformation [21j which corresponds to multiplying (or dividing because of 
the modular group symmetry r -s — > 1/r, i.e. the exchange of the two periods of the 
elliptic curve), the ratio of the two periods: 

k ^ kL= ^, r ^ 2r. (13) 

1 + fc 

The other transformationsyj correspond to r ^-)> A'' -r, for various integers A'^. However, 
in the natural variables of the model (as e^, tanh(ii'), k = s"^ = sinh {2K), but 
not the "transcendental variables" like r or the nome q), these transformations 
are algebraic transformations corresponding, in fact, to the fundamental modular 
curves. For instance, the Landen transformation (|13p corresponds to the genus zero 
fundamental modular curve 

j'-/' -(j+/)-(j-' + 1487.j-/+/2) 

+ 3 • 15^ • (16 f - 4027 jf + 16 j'^) (14) 

-12-SO^-iJ +j') +8-30^ = 0, 

which relates the two j-functions 

..,. ,,. (1 - fc^ + k^r .,, , ,. (i + l4p + fc4)3 

^^')= '''• (1-P)2.fc4 ' ^(^^)= ^^- (1-F)4.F ■ 

or to the fundamental modular curve: 

5%^^^ -l2-b^u^v^ • (u + w) +375 uw (16 u^ +16w^ -4027um) 
-M{v + u)-{v'^ + U%7vu +u^) +2i2-3^-uu= 0, (15) 
which relates the two Hauptmoduls u = 12'^/ j{k), v = 12'^/j(fci): 

t See for instance (2.18) in [2B]. 
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The fact that such transformations correspond to either multiplying, or 
dividing, t is associated with the reversibility^ of these exact representations of the 
renormalization group f25' . The "price to pay" is that these algebraic transformations 
are not one-to-one transformations, they are sometimes called "correspondences" by 
some authors. 

A simple rational parametrizatio njji of the genus zero modular curve (|15|) reads: 

1728 z 1728 2^ /212n 

" = "(") ^ (7TW " = (^T^^ = HvJ- (^^) 

Note that the previously mentioned reversibility is also associated with the fact 
that the modular curve (|15p is invariant by the permutation u ^-> v, which, 
within the previous rational parametrization (|16p. correspondfja to the Atkin-Lehner 
involution [30] z O 2^^/z. 

It has also been underlined in [201 HI] that seeing (IT51) as a transformation 
on complex variables (instead of real variables) provides, beyond k — 0, 1 (the 
infinite temperature fixed point and the critical temperature fixed point), two other 
complex fixed points which actually correspond to complex multiplication for the 
elliptic curve, and are, actually, fundamental new singularitietlSJ discovered on the x 
linear ODE [8l[33l[34]. In general, within the theory of elliptic curves, this underlines 
the interpretation of the (generators of) the renormalization group as isogenics of 
elliptic curves |25j . Hauptmodul^ modular curves and modular forms play here a 
fundamental role. 

Along this modular form line, let us consider the second order linear differential 
operator 

(z2 + 56z+1024) 240 



(z + 16)(z + 64) ^ z • (z + 16)2 (z -f- 64) ' 
which has the (modular form) solution: 

2^i([^^ J7]^M;1728 ,^ ' \ (17) 



12' 12^' ' " (z-f-16)3. 

256\-i/4 ^ / 1 5 ^ ^2 



-y ■.f.O^.y^l.M;!™ 



z + 16 / V42' 12' ' " (z + 256)= 

I The renormalization group (RG) is introduced, in the textbooks, as a semi-group, the 
transformations of the RG being non-invertible. We do not try here to address the embarrassing 
question of providing a mathematically well-defined definition of the RG. It is known that there are 
many problems with the RG, for instance for the (two-parameter) non-zero-temperature Ising model 
with a magnetic field [27! 1281, !29j. We just remark here, for a Yang-Baxter integrable model like the 
Ising model without magnetic field, that the RG generators, seen as transformations in the parameter 
space, identify with the isogenics of elliptic curves and are thus reversible in the r-space. 
ft Corresponding to Atkin-Lehner polynomials and Weber's functions. 

§ Conversely, and more precisely, writing 1728 z^ /{z + 256)'^ = 1728 z' /{z' + 16)^ gives the Fricke- 
Atkin-Lehner |30II31I involution z ■ z' = 2^^ , together with the quadratic relation z — z z' — 48 z z' — 
4096 z' = 0. 

^ Suggesting an understanding |21l 1321 of the quite rich structure |21| of the infinite number of 
singularities of the X ^ in the complex plane, from a Hauptmodul approach I21II32| . Furthermore, 
the notion of Heegner numbers is closely linked to the isogenies mentioned here |21| . An exact value 
of the j-function j{t) corresponding to one of the first Heegner number is, for instance, j = 12'^. 
f It should be recalled that the mirror symmetry found with Calabi- Yau manifolds 1351 1361 I571I38II39I 
I40| can be seen as higher order generalizations of Hauptmoduls. We can, thus, expect generalizations 
of this identification of the renormalization and modular structure when one is not restricted to 
elliptic curves anymore. 
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In fact the two pull-hacks in the arguments of the sam,e hypergeometric function are 
actually related by the fundamental modular curve (fT51) (see (HH)). Do note that, 
generically, the existence of several pull-backs for a hypergeometric function is a quite 
rare situation. The covariance (J17p is thus the very expression of a modular form 
property with respect to transformation r <-> 2 r, corresponding to the modular curve 

This example (see ([T7| ) is a simple illustration of the special role played by 
selected hypergeometric functions having several possible pull-backs (in fact an infinite 
number). This phenomenon is linked to elliptic curves in a deep and fundamental 
sense, namely the occurrence of modular forms and of isogenics represented by 
algebraic transformations, called by some authors "correspondence" , between the pull- 
backs. These algebraic transformations correspond to modular curves (and to exact 
algebraic representations of the generators of the renormalization group |25|). We will 
say in short, that such hypergeometric functions are "associated with elliptic curves". 
Simpler examples of isogenics associated with rational transformations, instead of 
"correspondence" like (IT5|) . are displayed in [25] . 



3. Modular form solution of F2 and the corresponding fundamental 
modular curve Xo(2) 

The exact expressions of the selected linear differential operators Z2, F2, and F3 
which emerged [5] in x^^\ can be found in |41], and the exact expressions of the 
selected linear differential operators L3, and L4 which emerged [7] in x^^\ can be 
found in |^. 

In the following we will not detail how we have been able to get the solutions 
of the linear differential operators F2, and F3, and the solutions of Z3, and L4. 
These details will be displayed in forthcoming publications. These (slightly involved) 
solution^a are displayed in .43^. We focus, here, on the structures, and mathematical 
meaning, associated with these solutions. 

Actually, the series expansion of the solutions of (the globally nilpotent) operator 
F2 gives more than G-series [44l|45]: it yields series with mte(;er coefficients, suggesting 
that F2 could also have, like the previous Z2, a modular form interpretation. 

The solutions of the second order linear differential operator |5] F2 can, actually, 
be written in terms of hypergeometric functions (the pi (x) 's are two rational functions) 

Piix)'/'- 2Fi([i, i],[i]; p{x)) + p2{xYl'- 2^i([^, J],[^]; p(x)), 
with the pull-back: 



1 (1 -4a;) (l-h6a;-hl3a;2 -^4x3) 



2 



p{x) ^ -—"- '-^ ^ '—. (18) 

^ ' 64 (l + 2x)'-x3 ^ ' 

A well-known symmetry of many hypergeometric functions amounts to changing the 
pull-back p(x) intqj q{x) = 1 — p(x): 

§ The reader can view all these linear differential equations, as well as their explicit 2^1 -type 

solutions, in |43j . 

X In addition to yix) — >■ q(x) = l/p{x). 
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where the selecteqjj singularities l + 3a; + 4 2;^ = 0, seen in x^'^\ and more 
specifically in ^2, clearly occur. 

Alternatively, the solutions of F2 can also be written as (the i?i's are rational 
functions) 

with two different possible pull-backs pi{x) and ^2(2^) reading respectively: 

{1 +8x +Ux'^ -36x^- 151 x^ - 188 x^ - 16 x^ - 64 x"^)^ 
1728-(l + 2a;)6. (1 + Ax)^ ■ {I - x)^ ■ {I + 3x + Ax"^) ■ x^ ' 
(1 + 8 a; + 14 a;2 - 276 x^ - 1591 x" - 3068 x"" - 1936 a;^ - 64 x'^)^ 
1728- (l + 3x + 4x2)2 . (l-a;)6. (1 +4a;)4. (1 +2a;)3-x3 ' 
Note that these two pull-backs are not related by a simple Atkin-Lehner involution |30] 
{p2{x) 7^ pi{N/x), with iV some integer). However, introducing the rational 
expression 

m) = -^—fp^^ (20) 

one actually finds that 

Pi{x) = ^(^)> P2(x) = R(q{x)). (21) 

The relation between these two pull-backs corresponds to the (genus zero) curve 
110592 • a2 /32 - 64 • [a^ + 13^) - 95232 ■ [a^ P + a P^) 
+ 6000 • (a2 + /32) - 1510125 • a jS 

- 187500 • (a +/?)+ 1953125 = 0, (22) 

with the simple rational parametrization deduced from (|2ip: (^(z), /3(z)) = 
{R{U), R{l/U)). One immediately finds that (HU is nothingbut ^ where (m, w) 
have been changed into (1/m, ^/v). Actually, up to a rescalin^ of f/ by a factor —64, 
the parametrization of the rational curve (|22p can be rewritten in a form where a 
o '2P' I z Atkin-Lehner involution is made explicit: 

i ( \ M w f ^ (-g + 16)' 1 (£+256)3n 

(«(.), /3W) - (17^^^' 1728^^^]' ^'') 

/3(z) = a(— ). (24) 

One recognizes (up to a 1728 normalization factor) the fundamental modular curve 
^o(2) 

^2^2 -{A +B)-{A^ + U87AB + B^) (25) 

- 40773375 ■AB+ 162000 • {A^ + B^) 

- 8748000000 ■ {A + B) + 157464000000000 = 0, 
and its well-known rational parametrization: 

.4.A,«.(?5^, B^a{^). (26) 

t Complex fixed points of the Landen transformation, Heegner numbers with complex multiplication 

of the elliptic curve |21l . 

^ R{U) = J2{-64:U)/1728 with J2{t) = (t + 256)='/1728/t2. 
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3.1. Dedekind eta function parametrization 

It is well known that the Dedekind eta function [461 I47j . in power 24, is a cusp 
automorphic form of weight 12, related to the discriminant of the elliptic curve. 
Recalling the Weierstrass' modular discriminant i48l, defined as 

A(r) = (2^)12 . r;(r)24, (27) 

which is this modular form of weight 12, we get rid of this (27r)^^ factor and define 

OQ 

A(g) - q-X{{l-q-f\ (28) 

n=l 

where q is the nome of the elliptic curve, that is, the exponential of the ratio of the 
two periodfia of the elliptic curve: q = exp(r). 

One can now introduce a "second layer" of parametrization writing the j-function 
as a ratio of Dedekind eta function 

\3 



n = 32{q) = ^^, A{n) = ^ -^^^. (29) 



A(g) (256 + J2) 

One deduces the alternative parametrization of ([25 



A 



(A, B) = (AiMq)), AiMq'j) = {a{^), ^(||^))> (30) 

making crystal clear that the fundamental modular curve (I25[) is a representation of 
T — )■ 2 • r, or g — )• g^ (and in the same timcp t -^ t/2). 

The Atkin-Lehner involutive transformation J2 ~^ 2^^/j2 and transformation 
q ^ q^ are actually compatible thanks to the remarkable "Ramanujan-like" functional 
identity on Dedekind rj functions 

4096 ■ A{q) ■ A{qy -A(g2)3 (31) 

+ {A{q) + 48 ■ A{q')) ■ A{q) ■ A{q^) = 0, 
yielding: 

(^) ^ -(^"/(|^))' '-) 

making ([^5|) and ([50)1 compatible. 

There are many other nice functional and differential relations on the Dedekind 
eta functions that are shown below. 

• The modular functions (see (19) in [49] page 16) 

^ / ??(6r)77(T) \i2 ^ 77(6r)S ??(t)^ 

V77(2r)7/(3T)/ ' ^ 77(2r)8 ?7(3t)4 ' ^ ^ 

have the following relation: 

« ^ »-ij^^ (34) 

This is exactly the covering necessary to see Z2 as a modular form (see equation (A. 3) 
ini). 

§ In the literature the ratio of the two periods of the elliptic curve often encapsulates a, 2 in (or i n) 
factor and one defines: q = exp(2i7rr). For some Calabi- Yau reason (see H125|l below) we prefer to 
write q = exp(T). 

f Thanks to the t •<-> 1/r symmetry of the modular group corresponding to the exchange the two 
periods of the elliptic curve. 
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• Differential equations are actually satisfied by modular forms [121 [5D] . 
Introducing the same t as in (155)) and the following function F{t) 

^ / r;(rM6T) y2 ^ (r?(2T)77(3r))^ 

Vr;(2r)r7(3T)/ ' {v{T)v{Qr)f ' ^ ^ 

one has an Apery's third order ODE ^ on the modular fornyl F{t). This ODE 
corresponds to the linear differential operatoitjj: 

(f2 - ut + 1) • t^ . Df + (6t2 - 153t + 3) • t • D'^ (36) 

+ {7t^ -112t +l)-Dt +(i-5), 

that reads in terms of the homogeneous derivative — t ■ d/dt: 

{f -Sit +l)-e^ + {3t^ ~51t)-e^ (37) 

+ {3t^ -27t)-0 +(<2 -5i), 

this operatoiLJj being linked to the modularity of the algebraic variety 



1111 




— \- y -\ t-zH \- w -\ = 


-- u, 


X y z w 





that is, to the one-parameter family of K3-surfaceo L53_: 

1 - (1 - XY) ■ Z - z-XY Z ■ {I- X)- {!-¥)■ {1- Z) = 0. 

• Another example of linear differential equations, satisfied by modular forms, 
can be found in page 18 of [i9] : 

V r][T)r](6T) / \ r]{2T)ri[bT) / 

the third order ODE on F{t) corresponding to the linear differential operator 

(64^2 +20 i +1)- 6*3 + {192t^ +30t)-e^ (39) 

+ {192t^ +18t)-0 +(64^2 +4i), 

• A third example [JU] is: 



/ 7y(3r)?](6T) x4 ^ (77(T)r?(2r))3 

Vr7fT)77f2r) / ' n(3T)v(6T) ' ^ ^ 



V r]{T)T]{2T) J ' 77(3r)77(6T) 

with a third order ODE on the modular form F{t). This ODE corresponds to the 
linear differential operator 

(81 1^ +Ut +l)-6^ + (243 1^ + 21 i) • 6^ (41) 

+ {243t^ +13t)-e +{8lt^ +3t), 

X F{t) and t are modular forms on ro(6) which has four inequivalent cusps oo, 0, 1/2, 1/3. F{t) 

is a weight 2 modular form. 

II In Apery's proof of the irrationality of f (3) a crucial role is played by the linear differential operator 

f Introducing the inhomogeneous order-two ODE corresponding to II37I I with the very simple rhs 
6 1, and considering the ratio of solution of I I37I I and of this inhomogeneous order-two ODE, one can 
build |49l a modular form of weight 4, by performing the third order derivative with respect to r, the 
ratio of two solutions of II36I I . 

^ The simplest example of Calabi-Yau manifolds are K3 surfaces. This Apery operator I I37II was seen 
in | 51| (see also [9|) to correspond to a symmetric square of a second order operator associated to a 
modular form (see also section 12.111 ). Along this line it is worth recalling that some one-parameter 
families of K3 surfaces can be obtained from the square of families of elliptic curves (see the so-called 
Shioda- Inose structures and their Picard-Fuchs differential equations ^52^ and see also relations (1.9) 
in [40]). 
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If one chooses two linearly independent solutions (Fi, F2) of these last order-three 
linear differential operators appropriately [15] , then i is a modular form of r the ratio 
of these two solutions. The function tir) satisfies a well-known third order non-linear 
ODE known as the Schwarzian equation in the litteraturqj: 

2Q{t)-[—) + {t, r} - 0, (42) 

where {z, t} denotes the Schwarzian derivative with respect to r: 

d d , , z(3) 3 // 2 

and where Q{t) is a rational function that can be simply deduced [19] from the 
coefficients of the fc-th order (here fc = 3) linear ODE on F{t). 

4. Modular form solution of F3, related to /ig, Apery and Z2 

The order-three linear differential operator F3 , occurring as a factor of the differential 
operator annihilating x , '^^is rationaly reconstructed in [S]. It can be seen to be 
homomorphic to the symmetric square of a second order operator. Similarly to what 
we had for F2, the (analytical at x = 0) solution of F3 corresponds to a series with 
integer coefficients, suggesting, again, a modular form interpretation. 

Actually the solutions of this second order operator can be expressed in terms 
of quadratic expressions of Legendre or 2F1 hypergeometric functions with a rational 
pull-back. The three solutions of F3 can be expressed in terms of Legendre functions 
where the pull-back Pi, in these Legendre functions, reads: 

= J_ (l^M(l±2^0(i±32xT (44) 

108 x'^ 

108 x^ 

The solutions can also be expressed as quadratic expressions of hypergeometric 
functions: 

Do note that the pull-back ((44]) is not unique. Another (Atkin-Lehner involution 
related) pull-back works equally well: 

Note that these two rational pull-backs are functions of x^ . These two pull-backs can 
be seen as a rational parametrization (a, b) = {Pi{x), P2{x)) of the (a, 6)-symmetric 
genus zero curve: 

- 625 + 525 • (a + &) + 3 6a + 96 • (a^ + b^) 

-528-{ba^ + b^a) +4-{a^ -\-b^) +4:32-a^b^ = 0. (49) 



J The simplest example of Schwarzian equation, associated with the complete elliptic integral K, 
corresponds to a rhs in II42I I reading — i ■ J^ i-TIii^ ■ '^^^ study of the Schwarzian equation is, in 
general, complicated, but Halphen found a simpler equivalent system of differential equations |54II55| . 
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Keeping in minqjj (|l5l) and (|48l) . we could have considered the algebraic curve 
relating {A, B) = (1 — Pi (x), 1 — P2{x)), which reads: 

-432^252 +4.(yl3^^3) + 336 • (A^ B + BM) (50) 

+ 381- AS - 12- (A^+B^) +12-(A + B)-4 = 0, 

which is rationally parametrized as [A, B) = {A{z), B{z)), where A{z) and B{z) 
read respectively: 

1 {z + Wf 1 (z + 64)3 

^(") - 1728^^' ^(") = 432^^' ^^1) 

where ^(2) and -B(z) are related by a Atkin-Lehner involution B{z) = A{2^'^ j z). 
This rational parametrization is extremely similar to the parametrization (j23p . (|26p 
of the fundamental modular curve (j25l) . One deduces, from (j5ip . the rational 
parametrization for the curve (H^ 

(z + 64)(z-8)^ (z + 16)(z-128)^ 

"^ 1728 ■ z ' ^= 432-^^ ■ ('2) 

where, again, 6(z) = a{2^^/z). Within parametrization (|52p . expressions 
(Pi(a;), ^2(2;)) (see (gll), (gT])) correspond to: z = -256 • x^ or z = -4/2;^. 

From ([5T|) it is thus tempting to interpret the new genus zero algebraic curve 
(|49p . or ((50)) . as a modular curv^ relating two Hauptmoduls corresponding to the two 
pull-backs (|44|) and (|47l) . similarly to what was found (see subsection (12. 1|) ) for the 
second order operator Z2 and its weight-1 modular form solutions. It was seen to be 
related [9] to a second order linear differential operator occurring in Apery's analysis 
of C(3): 

Ax- {x^ -Mx+l) ■D.J 

+ A-{1 -5lx +2x'^)-D^ +a;-10. (53) 



From its rational parametrization this new curve (|50|) is extremely similar to (|25p. the 
fundamental modular curve Xo(2). One has, of course, the well-known (and slightly 
tautological) algebraic geometry statement that all the genus zero curves of the plane 
are birationaly equivalent. But refereing to the "second layer" of parametrization (see 
pop above) can we say that this new curve is also a representation of r i — > N r and 
thus "truly" a modular curve ? 

4.1. The new curve i5(^) and the modular group ro(6) 

Seeking for hypergeometric functions with pull-backs that are not rational functions 
anymore, but algebraic extensions, we actually found another description of the 
solutions. The second order operator ()53p can be solved in terms of hypergeometric 

f And keeping in mind the well-known symmetry of many hypergeometric functions changing the 
pull-back p(x) into (((a;) = l—p(x). The other well-known symmetry Pi (z) •<-> 1/Pi(x) corresponds 
to x^ i — s. (1 -4a;2)/(l +32a;2)/4. 

X Modular curves of genus 0, which are quite rare, turned out to be of major importance in relation 
with the monstrous moonshine conjectures [561157) . In general, a modular function field is a function 
field of a modular curve (or, occasionally, of some other moduli space that turns out to be an 
irreducible variety). Genus means that such a function field has a single transcendental function as 
generator: for example the j-function. The traditional name for such a generator, which is unique up 
to a Mobius transformation and can be appropriately normalized, is a Hauptmodul (main or principal 
modular function). 
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functions 2^1 ( [I, 1], [§]; P± ) with the two possible algebraic (Galois conjugate) pull 



^L3' 3J' L2 

backs: 

1 a;4 - 23 a;3 - 156 x^ - 23 a; + 1 
^^ - -216 ^^ (^^) 



1 (.T-l)(a:^-7x + l) /--— 
^216 ^^ ^'-^^^ 



To get rid of the square root in (|54p. we introduce a parametrization of the rational 
curve y^ = 1 — 34 a; + a;^, namely: 

1 (w + 18)(u+16) 1288-u2 

2 u 2 u 

In terms of this rational parametrization (1551) the two possible pull-backs (|54l) read 
respectively: 

P,iu)^ _^(. + 24f(.^ + 12. 72)- 
+ ^ ^ 432 m-(u + 16)(m + 18)2 ' ^ ' 



1 (u + 12)2 (m2_ 48 u- 1152)2 ^ 288\ 

216 u2 (y _j_ iQ^2 (y ^ 18) - +V~y' 



or, more simply on ((5+(w), Q_(w)) = (1 — P+(w), 1 — P_(m)): 
„ , , _ 1 (^ + 12)*^ 

^+^"^ - 432^.(u + 16)(u+18)2' ^^^> 

n ( \ = n ('^\ - 1 (^ + 24)^ 

(^_^u; ^^+1 u J 216u2.(m + 16)2(u + 18)' 

Noticeably, these two pull-backs (|56p can actually be seen as another rational 
parametrization (a, 6) = (P+(m), P_(?i)) of the genus zero curve (|49p. The two pull- 
backs ([ST]) are another rational parametrization of the "new" algebraic curve ([5D)) . 
with (A, B) = (Q+(w), Q-(m)). 

Recalling the fact that the two pull-backs Pi{x) and P2{x) (see (|44|) . (|47|) ) were 
functions of x2, the correspondence between -Pi (a;) and P2(a;), and the two pull-backs 
([56]). corresponds to the following change of variables: 
_^iu + l6)-u^ 

128 u + 18 ^ ' 

The results ([5^ could have been obtained, alternatively, recalling the change of 
variable p4p . namely x = v • (1 — 9 w)/(l — u), which transforms the linear differential 
operator (1551) into (after rescaling by (1 — v)^^'^): 

2 I l-20z;+27z;2 _ l-3t; 

^ ^ (l-9^)(l-i;)i; " (l-9«)(l-«)z;' ^ ^ 

which corresponds to the (staircase-polygon |17| ) second order operator Z3 seen in 
Appendix A of [9] (see equation after (A.4)). The variable u is related to the previous 
variable u by u = 18 • (v — 1). 

However, recalling the weight-1 modular form interpretation ojj operator Z2 in [9] , 
we had the occurrence of Hauptmoduls Aiz = ^ (resp. ^4-) corresponding to the 
(genus zero [TB]) modular curvqjj 

$6(j,/) - Mf,j) = 0, (60) 

t Or of the order- two operator Z3 corresponding to staircase polygons 1171 . 

t Which amounts to multiplying, or dividing, the ratio of the two periods of the elliptic curve by 6. 
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obtained from the elimination of z betweeiijll: 

(z + 6)3(z3 + i8z2 + 84z + 24)3 

' = ^«(^) = z-{z + 9)Hz + 8r ' ^''^ 

z-(z + 8)3\ (z + 16)3 z-(z + 8)3 



■\ z + 9 J 



•''^Vz + 9/ z ° z + 9' 

. /z-(z + 9)2n (z + 27)(z + 3)3 z-(z + 9)2 

= -^s -T, — = ° -5 — > and: 

\ z + 8 / z z + 8 

' 2^ • 32 \ (15552 + 3888 z + 252 z^ + z^)^ (z + 12)= 

J6| 



J3 



z / z6(z + 8)2(z + 9)3 

, /z'''(z + 8)\ _ (z + 256)3 z3-(z + i 

(z + 9)3 / ^ ^ ° (z + 9)3 

, /z2.(z + 9)\ (z + 27) (z + 243)3 z^ • (z + 9) 



JM .;, 0,3 ) = .2 ° ., , nV3 (62) 



V fz + 8)2 y 



o 



(z + 8)2 J z3 (z + 8)2 



with the covering [9] 



72 a; 
" = (l-x)(l-4x) ' ^^^) 

which is a shght modification of ([8]). Let us introduce the alternative covering 
u = 2z, the two Hauptmoduls A4z = ■^ (resp. ^) read respectively: 

(6) ^ 110592 ■ z. ■ (u + 16)3 (u + 18)^ 

1 ^ ^ (12 + u)3-(192 + 336u + 36m2 + u3)3 ^ ' 



^o3/^ (" + 32)% ^ ^ ^.(^ + 16)3 ^ 
V^^/V 4u JJ U.(u + 18)J' 



Pf(.)= Pr(^) (65) 



{u + 18) 

3456 -u*^ • (u + 16)2(ii + 18)= 
{u + 24)3 (124416 + 15552 u + 504 u^ + u^f 

fl23/f^" + ^^^)'^^^ "'('' + ^^^ 






2^2 yy (w + i8)3 

The relation between the Apery operator ([5^ and the Z2 weight- 1 modular forms 
associated with (|B(7|) . seems to say that there should be some (at first sight totally 
unexpected) relation between hypergeometric functions with the pull-backs (j56p and 
the hypergeometric functions with the pull-backs (|M)) . We have actually been able to 
find such a "quite non-trivial" relation 

= 2i/2.p.C+(z.). 2Fi([| ^],[^]; P+in)) (66) 

-p.C_(^).2Pi([|^],[^]; P-{u)), 
where C+{u), C-{u), C6{u) and p read respectively: 

/(U + 24)2x1/2 / 16^2 x2/3 /u2_|_;^2u-72 

V u / ' \{u + 16) (u + 18)2/ ' V 64^^^ 

^ Here " o " denotes the composition of two rational functions i.e. f{z) o g(z) = f(g{z)) 
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'2- (w + 12)2x1/2 / 4u \2/3 /w2_48u_ii52x 

u ) \{u + 18) {u + 16)2/ 'V 16u /' 

144 u2 Ni/4 2 (r(2/3))^ 



Vl82- (12 + u) (192 + 336 u + 36^2 + ^3)7 ' 3 

Note however, that Pi (u) (or Pj (")) cannot &e ea;j)resse(i as a rational function 

of the two pull-backs P±{u) (see ^^). The relation between ae = P{ {u) and 

a = P+iu) (resp. fee = P2 (u) and 6 = P-{u)) reads a (necessarily genus zero) 
algebraic curve 

{a~l)-{9a-25f -al (67) 

+ 8 (a- 1) -(1458 02 -1215 a +125)- ae + 16 = 0, 
and the same genus zero algebraic curve where one replaces {ag, a) — >■ {be, b). Note 
that the relation between P| '{u) and P-{u) (resp. Pj (^) ^^'^ P+i'^)) i^ niuch more 
involved. 

Recalling (|60p one deduces, from the previous calculations, that the (genus zero) 
modular curve 

-123 12^, 



\ ae be / 

is actually "equivalent" to the genus zero algebraic curve (|49p up to the algebraic 
covering (|67|) . 



.^.^. Dedekind parametrization of the new curve 

Let us revisit the rational curve ([50l) that we rewrite (with a rescaling of A and P by 

1728): 

- 2/2^2 _^ 16 (j^ ^ ^) (^2 + 83 2/z + y2) _ 82944 • (.z2 _,_ y2^ ^gg^ 

+ 2633472 • yz + 143327232 • (y + z) - 82556485632 = 0. 
The curve (j69p is rationally parametrized by: 

J2 ^J2/ 16 -72 

The correspondence with the previous rational parametrization is j2 = 4-z. Similarly 
to what was done for the fundamental modular curve (j25[) . we can introduce a second 
"layer" of parametrization, writing J2 as a ratio of Dedekind eta function ([28]) 

■ ( \ A ^(«) 4r ^ (256 + J2)' 
^2 = J2(<Z) - 4^^, A(j2) = —-2 , 

which yields the following parametrization for (|69p (A(j2) is the same function as in 

which makes clear that the algebraic curve (l69l) is a representation of g — > q^. The 
compatibility between the Atkin-Lehner involution J2 ^^^ "^^^ I H and the q ^ q^ 
transformation, corresponds to 

<W^ -' ^(^"/(4^)) (70) 
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which is nothing but ([32)) corresponding to the functional equation (|3T|) . 

Matching ?;/1728 or z/1728 (see dH]), gT])) with 1 - Pi or 1 - P2, one gets: 

j2 = -1024x2 or: j^ = -^. (71) 

yielding a straight interpretation of the x variable in the 71-fold integrals of the Ising 
model in terms of Dedekind eta function, and more precisely, of the discriminant A: 

5. Modular form solution of L3: is it associated with the new curve (|50p 
or with the fundamental modular curve Xo(2) ? 

In [7] an order-three linear differential operator L3 was found as a factor of the minimal 
order operator for x*'^^- This order-three operator is (as it should) globally nilpotent, 
and one can see that it is reducible to an order-two operator in the sense that it is 
homomorphic to the symmetric square of an order two linear differential operator: 

x-{l -16xf-{l -Axf -Dl (73) 

+ (1 - 24 x) (1 - 16 a;) (1 - 4 xf ■ D^ + 108 a;^ 

yielding solutions in terms of hypergeometric functions like 2-Fi([l/8, 3/8], [1/2]; Vi{x)), 
with the pull-back Vi{x) reading 

^^(^) = (1-1mT-^'4.)2 ' (^4) 



or more simply: 



l-V,{x) ^ - f' f (75) 

(1 — 162:) (1 — 4x)^ 



5.1. An "uneducated" guess 

A simple calculation shows that, some "Atkin-Lchner" transform of ([71]) 
^ . X ^ / 1 \ 4- (3- 16x)2 -a; 

or more simply 

1 - ^2^ = n T^^ TTT^' (^^) 

(1 — 4a;) (1 — 16a;)^ 

provides another rational parametrization (a, b) — (Pi (a;), 7^2(2;)) of the new genus 
zero curve (|49p . 

Let us introduce the x^-dependent transformation 

one actually finds a nice relation between the two pull-backs Vi (x) and this "guessed 
candidate" 7^2(2;), for the other pull-back of L3 (if any ...): 

Pi(i?(x)) = Fi(x), V2[r{x)) = P2{x). (79) 
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These two equalities are actually compatible with the two Atkin-Lehner involutions 
for {Vi, V2) and (Pi, ^2)1 because of the nice functional relation on R{x): 

r(—^ = ^-^^. (80) 

\8xJ 64-i?(x) ^ ' 

It is, thus, extremely tempting to imagine that L3 is, like F^, related to the new 
modular curve P^j) or ((5(I)) . Furthermore, this would yield some (deep ...) relation 
between the singularities of the x^^"'^ and singularities of the ;^(^"+i), that is to say, 
between the low and high temperature singularities of the susceptibility of the Ising 
model. This is not the case: L3 is, in fact, related to the fundamental modular curve 
^o(2). 

5.2. Modular form solution of L3: the fundamental modular curve Xq{2) 

Actually, leaving, again, the framework of rational pull-backs, one gets the two (Galois- 
conjugate) algebraic pull-backs 'P±[L3]: 

3456 • V±[L3] = (81) 

(40a;2 - 17.T + l) (400 x" - 928x=' -I- 297a;2 - 31 a; -Fl) 



X' 



6 



(1 -12a;)(l -4x)(l-7x)(25a;2-i7a;+l 



± ^^ — — t: — ■ \/l - 16x. 

x" 

The relation between these two Galois-conjugate pull-backs actually corresponds to 

(j22|) which is nothing (up to a 1728 rescaling factor, see ()25|) ) but the fundamental 

modular curve Xo(2). 

To get rid of the square root singularity, we introduce the variable y: 

y2 = 1 -16x, i.e. ^ = -T^ (2/'-l)- (82) 

16 

The two previous algebraic pull-backs become respectively: 

J_ {5y^-9y^ + 15y-3f J_ {5y^ + 9y^ + 15y + sf 

27 (j/ + 1)6 (y - 1)3 ' 27 {y + lfiy-lf ' ^ ' 

which can be seen to be a rational parametrization of (I22p . Recalling the previous 
parametrization (j23p one finds that the z variable in (j23p . must be equal to z = 
64 ■{y + lf/{y - if, or z = 64 • (y - lf/{y -t- 1)^. 

6. The puzzling L4: preliminary results on £4 

Let us now focus on the order- four linear differential operator L4, discovered as 
a factor of xlf^ and that we were fortunate enough to get exactly by rational 
reconstructiorlj [?]■ Let us display a few results on L4. 

6.1. Negative results on L4 

Suppose that L4 is equivalent (in the sense of equivalence of linear differential 
operators [TOllTI]) to a symmetric cube of a second order linear differential operator 
L2. Take a point x — a, and suppose that the highest exponent of ln(a; — a) that 

I For details on the rational reconstruction see | 58| . 
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appears in the formal solutions of ^2, at x — a, equals p. Then the highest exponent 
of ln(a; — a) that appears in the formal solutions of L4, aX x ~ a, must be 3/9. Now 
look at the formal solutions of L4, at a; = 1/8, one gets a contradiction. Similar 
reasoning (using both x = 1/8 and x — Q) shows that L4 can also not he related to 
the symmetric product of two second order operators. 

Our preliminary calculations show that L4 is not 4i^3-solvable if one restricts to 
rational pull-backs, in the sense that there is no 4F3 differential operator that can 
be sent to L4 (with a change of variables x — > rational function in x, followed by 
multiplying by exp( J (rational function in x)), followed by Homomorphisms) to L4. 

In essence, if we allow the following functions in C(x): exp, log, and any pFq 
function, composed with any rational functions, and anything one can form from 
those functions using addition, multiplication, derivatives, indefinite integral, then we 
believe that L4 is not solvable in that class of functions. 

6.2. Positive results on Li 

The order-four linear differential operator L4 exhibits, however, a set of nice 
properties. Let us display some of these nice properties. 

• An irreducible linear differential equation (resp. irreducible linear differential 
operator) is said to be of Maximal Unipotent Monodromy (MUM) if all the indicial 
exponents at x = are zero (one Jordan block to be Maximal) . The formal solutions 
of an order-four MUM linear differential operator can be written as 



(84) 

ys- 

u z, 

The indicial exponents of L4 at x = 0, read —6, —4, —4, and 0. Therefore, the 
order-four operator L4 is not MUM, however the formal solutions can be cast exactly 
as for a MUM linear differential operator. The formal solutions of L^ can be written 
as 

yo = yo, 

yi = yo •ln(x) +yio, 

1/2= -yo •ln(a;)2 +^21 •ln(x) +2/20, (85) 

ys = ^ yo • ln(.T)^ + 2)32 • ln(x)^ + ysi • ln(x) + 7/30 


Such particular form for the set of series solutions is often obtained for irreducible 

operators. 

There are four independent series that can be chosen as yo, yio, 2/20 a-nd ^30. The 

other series in front of the log's should depend on these four chosen series. For L4, 

these series read 

2854486264697 

y'' = y'' ^ 459375-106 ^"' 

.1. 38100003421933 

y^' - 2^'" + 11484375-104 ^°' 



yo = 


yo, 


yi = 


yo -ln(x) +yi. 


y2 = 


-yo • ln(x)^ + yi - ln(a;) + j/2, 


ys = 


-yo - In(x)-'' + -yi - ln(x)^ + y2 ln(x 



yo ^ 


yo, 


yi = 


yo •ln(a;) +yio, 


y2 = 


-yo • ln(a;)2 + yw ■ \n{x) + 


ya = 


-yo • H^f + ^yio ■ ln(a;)2 
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. _ . 85327 , 103266884399422867 
ysi - y2o +^^yio 504 . 10" ^°' 

where we see that it is a matter of combination to cast the formal solutions (|85p in the 
form ([84]). Making the combination j/2 — C21 yi and ya — 2cz2y2 — (C312 — 2c2iC32)yi, 
(where C21, resp. C32 and C312, are the coefficients in front of j/o hi ^21, resp. ^32 
and y3i, in ([86])), one obtains the new set 



2/20 - C2iyio, (86) 



+ (y20 + C3iiyio - 2c32Cioyio) • ln(a;) 

+ 2/30 - 2c32 • 2/20 + (2C32C21 - C312) • yio, 

where (y20 -C2iyio) identifies with (^20 +C3iiyio - 2c32Cioyio) since C21 = -C311 + 

2C32C10 for the actual values of the combination coefficients. 

The remaining difference between a MUM linear differential operator and the 

formal solutions of L4 is the beginning of the series of the non leading log's. For 

instance the local exponent —4 being twice, one should have a series starting as x^^ 

in front of the log's like yio 

,1 11 9329 8023 99922803261913 

^^° ~ 84000 x^ ^ 16800x3 ^ 336000 x^ ^ 8400~r ^ 3675000000000 ^ " ' 

We may imagine that by acting by an intertwinner on the formal solutions, one ends 

up with series starting at x, i.e. L4 may be equivalent to a linear differential operator 

which is MUM. 

• The order- four linear differential operator L4 is of course globally nilpotent [5|. 
The p-curvature [59] of the order-four globally nilpotent differential operator L4 can 
be put in a remarkably simple Jordan form: 

10 

10 

1 





(87) 



Its characteristic and minimal polynomial is T'*. Such an operator cannot be a 
symmetric cube of a second order operator in C{x)[Dx]. We, however, determined the 
differential Galois group 60, 61 of this linear differential operator L4 and actually 
found that it ij^jj the symplectic group SP{4, C). A crucial step to exhibit this 
symplectic structure, amounts to calculating the exterior square of L4, and verify that 
this, at first sight order-six, exterior square either reduces to an order-five operator, 
or is a direct sum of an order-five operator and an order-one operator with a rational 
solution. Li corresponds to this last scenario. 

ft SP{A, C) is not the monodromy group: it is equal to the Zariski closure of the (countable) 
monodromy group, i.e. the differential Galois group. 
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Denoting yi the four solutions of an order-four linear differential operator, the 
exterior square of that operator is a linear differential operator that annihilates the 
expressions 

w^,J = y^ 5 - % ^' * ^ .?■ = 0' 1' 2, 3. (88) 

It should, then, be (at first sight) of order six. 

The six solutions ([55)) of the exterior square of a MUM order-four operator, 
contain log's with degrees (at most) 1, 2, 3, 3, 4 and 5. There are then two solutions 
(with the same degree in the log's) that can be equal 

d d d d 

yo -i-ys - ys -rvo = yi -r-v^ ~ y^ -rv^- (89) 

dx dx dx ax 

When this happens the exterior square, annihilating five independent solutions, 
will be of order five. This is equivalent (see proposition 2 in [HI]), for the coefficients 
of the order-four linear differential operator, to verify condition (J122I) below. 

Here, computing the exterior square of the linear differential operator £4, one 
finds an order-six linear differential operator with the direct sum decom,position 

ext(2)(L4) = Ni®N5, (90) 

with 

sol{N,) = ^, (91) 

D{x) 

N{x) = -12 + 2548 X - 502593 x^ + 43407720 x^ - 1959091320 x^ 
-^52738591890 x^ - 904049598675 x^ + 10126459925120 x^ 
-74115473257440 a;® + 350453101085400x^-1133589089074624x1° 
+4059589860750336 x" - 25595376023494656x^2 
+141123001405931520 x^^ - 440315308230574080 x" 
+705909942330064896x^5 - 496507256028790784x^6 
+140082179425173504 x", 

D{x) = x^ • (1-16 x)^'''- (1-4 x)2 (8-252 x + 1678x2 -3607x3 -4352 x"*). 

This decomposition induces on the six solutions ([88)) of the exterior square, the 
following relation (up to a constant in sol{Ni)) 

9701589902493^04 + 609600054750928 u;o,2 + 91875 • 10^ • (wi,2 - ^0,3) = sol{Ni), 

which shows the occurrence of a relation between the four solutions yo, yi, 2/2 and j/3 
of L4 and their first derivative. 

• Along this globally nilpotent line, the L4 operator is "more" than a G- 
operator [331|3S], with its associated G-series. The series solution (analytical at x = 0) 
sol{L/i) is a series with integer coefficients in the variable y — x/2: 

sol{Li) = 175 + 34398 y + 4017125 y^ + 362935156 y^ 
+ 28020752579 y^ + 1943802285620 y^ + 124761498220195 y^ 
+ 7549851868859190 y^ + 436341703365296321 y^ 
+ 24309515324321362986 y^+ 1314618756208478845353 yi° 
+ 69377289961823319909960 y" + 3588051829563766082490527 y^^ 
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+ 182471551181260556637299032 yi3 

+ 9150139649421210256395488775 y" 

+ 453470079520701103056020155546 y^^ (92) 

+ 22252827613097363700809754930653 y^^ 

+ 1083008337798028206538475233669454 y^'' 

+ 52344841647844780032111214432202429 y^^ 

+ 2515396349437801867561610046658122428 y^^ 

+ 120295197044047707889910797105191140059 y'^° 

+ 5729990034986443499765238359785037134524 y^^ 

+ 272033605883471055363216581302378024952171 y^^ 

+ 12879727903873470148364481804530391226578654 y23 + ... 

7. Various scenarios for L4 

The order-four linear differential operator L4 is thus slightly puzzling. It has clearly a 
lot of remarkable properties but cannot be reduced in a simple, or even in an involved 
way (up to rational pull-backs, up to operator equivalence [inilll] i-e. homomorphisms, 
and, up to symmetric powers or products), to elliptic curves or modular forms. Is this 
operator going to be a counter-example to our favourite "mantra" that the Ising model 
is nothing but the theory of elliptic curves and other modular forms ? 
Let us display a few possible scenario^. 

7.1. Hypergeometric functions, modular forms, mirror maps 

Because of the globally nilpotent character of L4, we have some "hypergeometric 
functions" prejudicqj. Furthermore, we would also like to see some "renormalization 
group" symmetries acting on the solutions, may be in a more involved way than 
what was described, in subsection (|2.2|) and displayed in [25] , as isogenics of elliptic 
curves. We seem to have obstructions with rational pull-backs on hypergeometric 
functions. We should therefore consider generalizations to algebraic pull-backs, but 
we expect the pull-backs to be "special" possibly corresponding to modular curves. 
Furthermore the integrality property of the solution series ()92p (integer coefficients for 
the series) suggests to remain close to concepts, and structures, like modular forms, 
theta functions (which are modular forms) and mirror maps j39[ I54j . 

Along these lines, let us recall a selected hypergeometric function closely linked 
with isogenics of elliptic curves, modular form and mirror maps. Let us recall (see 
comments after formula (6.2) in 'ET) that|jj 



3-F2 



([y4Ml.ll;.) (93) 

= (l-.)--,f,([i2,il,[Ml;jj^ 



^ Taking into account the previous known results on L4, and consequently having some overlap. 
X It has been conjectured by Dwork 63 that globally nilpotent second order operators are necessarily 
associated with hypergeometric functions. This conjecture was ruled out by Krammer 1641 1651 for 
some examples corresponding to periods of abelian surfaces over a Shimura curve. 
t We have studied in some details the renormalization transformation z — > —4z/{l — z)^ in 1251 . 



Ising model, modular forms and Calabi-Yau 22 

is a modular map 62 . More generally, one has remarkable quadratic relations (see 
(6.3) in [62] )■ where the Landen transformation 

' ^ ITT^ (94) 

(1 + zY 

clearly pops out. 

This selected hypergeometric function satisfies the quadratic relatior|§|: 

3F2([i, i, !],[!, l];4t.(l-t)) = (2^i([i,i],[l];t))'. (95) 

which yields relation ^49l|67] (g denotes the nome): 






which is naturally associated with a very simple mirror map involving theta functional 



'M ?M fl4, ^ g dz 



-('^) = 4 i^ ■ WT' (^til)= —^= ■ ^. (97) 



Generalizations of this kind of relations are certainely the kind of relations we are 
seeking for, for the solutions of L4, but, unfortunately, this requires a lot of "guessing" 
(of the selected hypergeometric function, of the mirror map, of some well-suited ratio 
of theta functions, ...). 

7.2. Hadamard product 

In our "negative" result section (j6.ip we saw that L4, which is a globally nilpotent 
operator, cannot be reduced to 2F1 hypergeometric functions associated with elliptic 
curves (modular forms), up to transformations compatible with the global nilpotence, 
namely equivalence (homomorphisms) of linear differential operators, symmetric 
powers, or symmetric products of linear differential operators (even up to rational 
pull-backs). Having in mind this idea of compatibility with the global nilpotence, one 
more operation, a "product" operation, can still be introduced, namely the Hadamard 
product which quite canonicall^ builds globally nilpotent differential operators from 
globally nilpotent elementary "bricks". 

It has been seen in G. Almkvist and W. Zudilin [62], that one can build many 
(Calabi-Yau) order-four operators from the order-two elliptic curves operators, using 
the Hadamard product of series expansions. The Hadamard product it|jj, just a 
convolution [69] (|z| < \w\ < 1): 

f*9iz) = ^- j f{w)-g{z/w)-—. (98) 

Let us recall the order-two (elliptic curve associated) operator 

Le - x-{l-x)-Dl + {l-x)-D, + i, (99) 

§ Compare this relation witii tlie Bailey theorem of products |66| : {iFiii-nt 2\>['^\'y^))'^ — 

4F3([i,i, i,l],[l, 1, l];4t.(l-i)). 

tj Which are, as we know, modular forms of half integer weight. 

§ The fact that the global nilpotence is preserved by the Hadamard product is a consequence of the 

stability of the notion of G-connection under higher direct images for smooth morphisms ,68|. 

t Deligne's formula (simple application of the residue formula). 
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which has EllipticE{x^''^) as solution, and let us consider the Hadam,ard product 
of (the series expansion, at x = 0, of) EllipticE{x^/'^) with itself. This 
Hadamard square of (the series of) EUipticE is actually (the series of) a selected 
4i<3 hypergeometric function 

2 • EUipticE 2 • EUipticE 



- 4^^ 



3([-l/2, 1/2, 1/2, -1/2], [1,1,1]; x), 



which is a solution of the globally nilpoten\^ fourth order linear operator that we will 
write Had^{LE)'- 

Had'^iLs) = -1 -8-(a;-2)-i:>^ + 8 • (14 - 13 a;) • x • Z?^ 

+ 96-{l-x)-x^ -Dl + m-{l-x)-x^ ■ D'^. (101) 

The Jordan form of the p-curvature [9] of the globally nilpotcnt fourth order 
linear operator (jlOip actually identifies with the 4x4 matrix (|87p. of characteristic 
and minimal polynomial T^. Such a linear differential operator cannot be a symmetric 
cube of a second order operator in C{x)[Dx]. We can, however, certainly say that 
the globally nilpotent fourth order linear differential operator (jlOll) is "associated 
with elliptic curves", and we will also say, by abuse of languagqj, that the linear 
differential operator (jlOip is the Hadamard product (at x = 0) of (|99l) with itself, 
or the Hadamard square of the linear differential operator (|99|) . and we will write 
Had^{LE) = Le * Le- Several examples of "Hadamard powers" of the complete 
elliptic integral K are given in ( [Appendix A.l[ ). 

In our miscellaneous analysis of various (large order globally nilpotent) linear 
differential operators, we try to decompose these (large) operators into products, 
and ideally direct-sums [IJ El H [70], of factors of smaller orders. We then try, in 
order to understand their "very nature" , to see if these irreducible factors are, up to 
equivalence of linear differential operators, and up to pull-backs, symmetric products 
of operatortljj of smaller orders. Since the Hadamard product quite naturally builds 
globally nilpotent operators from globally nilpotent ones, and since it already provided 
examples j71| of (Calabi- Yau) order-four operators for which the corresponding 
mirror symmetries are generalizations of Hauptmoduls (basically products of elliptic 
curves Hauptmoduls 150)). we can see the Hadamard product as a quite canonical 
transformation to add to the symmetric product of linear differential operatorSj^. 

We will say that an irreducible differential operator is "associated with an 
elliptic curve" if it can be shown to be equivalent, up to pull-backs, to a symmetric 
product, or a Hadamard product, of second order hypergeometric differential operators 
corresponding to elliptic curves (see PS]). If the differential operator is factorizable, we 
will say that it is " associated with an elliptic curve" , if each factor in the factorization 
is. 

tl The Hadamard product of two hypergeometric series is of course a hypergeometric series. The 

minimal operator of a G-series is globally nilpotent, and the Hadamard product of two G-series is a 

G-series. 

f The Maple command gfun[hadamardproduct](egi, 652) returns the ODE that annihilates the 

termwise product of two holonomic power series of ODEs, eqi and eq2. 

ft That is simple products of the solutions. 

X For operators, not necessarily irreducible, this amounts to considering five "Grothendieckian" 

operations: three products, the products of the operators, the products of the solutions of the 

operators (symmetric product), the Hadamard product (convolution, Fourier transform), as well 

as the operator equivalence, and the substitution (pull-back) |44l 1451 159) . 
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Is L4 in [7j an operator "associated with an elliptic curve"? This looks like a 
quite systematic (almost algorithmic) approach. In practice, it remains, unfortunately, 
(computionally) very difBcull|§| to recognize Hadamard products up to homomorphism,s 
transformations . 

7.3. Calabi-Yau and SP{A, C). Recalling three-fold Calabi-Yau manifolds 

We have discovered a symplectic 5'P(4, C) differential Galois group for L4. Many 
order-four operators (often obtained by Hadamard product of second order operators) 
and corresponding to Calabi-Yau ODEs, were found to exhibit a symplectic S'P(4, C) 
differential Galois group, to such a large extend that it may be tempting, for order-four 
operators, to see the occurrence of a SP{A, C) differential Galois group as a strong [75] 
indicatiorU in favour of a Calabi-Yau ODE (75]. On the other hand, one may have 
the prejudice that Calabi-Yau ODEs and manifolds, which are well-known in string 
theory, have no reason to occur in (integrable) lattice statistical mechanics. This is 
no longer true after Guttmann's papeilj which showed very clearly ,74^ the emergence 
of Calabi-Yau ODEs in lattice statistical mechanics. 

At this step, let us recall the famous non-triviaUl example [35j of Candelas et 
al. of (three-fold) Calabi- Yau manifold. The order-four linear differential operator (in 
terms of the homogeneous derivative 9 — z ■ d/dz) 

6|4 - 5 z • (5 6* -h 1) • (5 6* -f 2) • (5 6* -f 3) • (5 61 + 4), (102) 

has the simple hypergeometric solution: 

4F3([i,||^],[l, 1, l];55z), (103) 

which is associated with the three-fold Calabi-Yau manifold [35j [75] : 

x\ +xl +xl +x\ +xl - Z~^/^ • Xi X2 X3 0:4 2:5 = 0. (104) 

Actually the hypergeometric solution (I103|) can be written as a multiple integral with 
an algebraic integrand having (|104p as a divisor. Being a multiple integral with an 
algebraic integrand it is, in mathematical language [9 "a period^ , and, consequently, 
the associated order-four linear differential operator (jl02|) is necessarily [9] globally 
nilpotent. 

The differential Galois group of (I102p is actually the symplectic [7?, ^75] group 
5*^(4, C). More precisely, the Picard-Fuchs linear differential operator (|102p . with its 
solution (|103|) [x — z), reads: 

x^ ■ (l-3125a;)-i:';^ + 2a;3 • (3 - 12500x) • L>3 

+ x^ ■ {7 - AbQQQ x) ■ Dl -f x- (1- 15000a;)-i:»:r +120a;, 
which can be written in the form (3.9) in [35] . when rescaling z — h^ ■ x: 

d^F{z) (4z-3) d^F{z) 1 (72Z-35) d^F{z) 

dz'^ ~ (1 - z) • z ' rfz3 ^5 (1 - z) • z2 ' dz^ 
1 (24 z- 5) dF{z) 24 F{z) 

~ 5 (1 - z) • z3 ' dz ~ 625 (1-z) -z^' 

§ The transformations "Homomorphisms" and "Hadamard product" mess up each other quite 

badly so a Hadamard product becomes difficult to recognize after homomorphisms (i.e. gauge) 

transformations . 

tl In fact Calabi-Yau ODEs do not reduce to SP{A, C) differential Galois group. 

t See [73| for the Fast Track Communication. 

^ Elliptic curves can be seen as the simplest examples of Calabi-Yau manifolds. 
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Along this line, some list of Calabi-Yau ODEs and Calabi-Yau linear differential 
operators have been obtained [7T] by G. Almkvist et al. seeking systematically for 
order-four differential operators obtained from Hadamard product constructions of 
second order operators, often within a symplectic and MUM framework. Such long, 
and detailed, list of Calabi-Yau differential operators are precious, but, again, it is not 
straightforward to see if an order- four operator, like L4, reduces to one of the Calabi- 
Yau differential operators in such lists, up to homomorphisms, and up to pull-backs. 

7.4- The 4-^3 scenario 

As far as order-four operators that cannot be simply reduced to elliptic curves are 

concerned, we already saw J33J, in the Ising model, an example corresponding [TTj to 

the form factors C^^'{k, n), expressed in terms of a 4-F3 hypergeometric function: 

,,, , ^/A+k + n 1 + k + n 2 + k + n 2 + k + n^ 
b{k,n)^ 4F3([ , , , ], 

[l + k,l + n,l + k + n];16x). (105) 

It is solution of an order-four linear differential operator which can be written, in 
terms of the homogeneous derivative 9 (in the usual quasi-factorized form for nF-a-i 
hypergeometric function): 

J,„= 16...(.+i±A±^)^(^+^±i±^)^ (106) 

- [e + k) ■ {9 + n) ■ [9 ^ k + n) ■ e. 

All these operators (|106|) are, in fact, homomorphic (see ( [Appendix B[ )). The 
linear differential operator ()106|) is not MUM (except for k — n — Q), however 
b{k,n) is clearly a Hadamard product (see (G.l) in |33) ) of two algebraic functions for 
k and n integers (or an algebraic function and a 3-F2 function otherwise). 

The exterior square of Jfc_„ is a sixth order operator which is invariant by A; -f-)- n. 
Noticeably, this exterior square of Jfe_„ has a very simple rational solution: 

- where: P= (1 - 16 a;) • x''+"+\ (107) 

which shows that one actually has a symplectic structure when k + n is an integer 
number. Actually, performing the direct-sum factorizatiorQ of the exterior square of 
Jk,n, gives (when k 7^ ±n), with P given by (|107|) 

Ext^Jk^r.) = ^i^^ © ^1^^ ® (q'^^-Q'^^], with: (108) 

ax ax \ Fk,n 

where Pk^n is a polynomial, where N and Af are integers depending of k and n, and 
where Q2 and Q2 are equivalent, and homomorphic to Q2 {k = l,n = 0) 

^9 „ (3-64x) 3-98x , , 

(1 — 16 a;) • a; (1 ~ lo x) ' ^ 

which solutions can be expressed in terms of hypergeometric functions: 

;^-.^i([yU2];l-16 

a;2.(l-32a;)3/2- ^^^([4' 4^' W' (1^32^) ' ^^^^^ 

t DFactorLCLM in Maple. 
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The integer M is equal to if fc — n is even, and is equal to 1 if fc — n is odd, 
the integer N and the degree of the polynomial reading respectively: 

3, ,, 1, ,, _, n + k \n-\-k\ 3 1, _,,„_i. 

- -{n + k) -\ \n-k\ +1, h • (-1) • 

2 ' 2 ^ ' 2 2 2 2 ' 

The symplectic form of the exterior square of Jfe.„ is singular if, and only if, 
k = ztn. The exterior square i?a;i^(Jfc,„) has no direct sum factorization for k — ±n. 
It factorises in the product of an order-one, two order-two and an order-one operators. 

Furthermore, the function 

a{k,n) - (^^^^''yb{k,n), (111) 

which corresponds to the form factor C^^'(k, n), has a series expansion with integer 
coefficients: 

fk + n\ (k + n+l)(k + n + 2)'^ fk + n^ 
"(''") ^ ( k j+ (n + l)Wl) ■( k 

aoik, n) /k + n\ ^ a^iik, n) /k + n\ r, 

+ ^y^-( k )■' ^^-^i k j-^ +••• 

where 02 (fc, n) and 0:3 (fc, n) read respectively: 

(fc + n + 1) (fc + n + 2) (A: + n + 3)2 (fc + n + 4)2 
{k + l)(k + 2){n+l){n + 2) ' 

{k + n+l){k + n + 2){k + n + 3)(k + n + 4:'f {k + n + 5f {k + n + 6y 
(fc + 1) (fc + 2) {k + 3) (n + 1) (n + 2) (n + 3) 

8. i4 is Calabi- Yau 

8.1. Warm-up: discovering the proper algebraic extension for the pull-backs 

The 4^3 function satisfies a linear differential operator L4^3 with three singularities 
0,1,00. 

The singularities of L4 at x — 1/16, and at x = 00, have exponents: 
integer, integer, half-integer, half-integer, and have only one logarithm there. This 
configuration is not compatible with any of the singularities of £4^3 under rational 
puUbacks. 

The singularity at x = 1 of ^4^3 has exponents 0,1,2, A, where A depends on the 
parameters of the 4F3 function. The exponents 0, 1, 2 correspond to solutions without 
logarithms. Thus, by choosing the 4F3 parameters to set A to an integer (we take 
A = 1) we get one logarithm at x = 1. Then the x — 1 singularity of L4_3 has the 
same number of logarithms as the x = 1/16 and x = 00 singularities of L4. However, 
under rational puUbacks there is still no match because the exponents of £4^3 at x = 1, 
which are now 0, 1, 1, 2, do not match (modulo the integers) the exponents of L4 at 
X = 1/16 and x = 00. 

A pullback x i— > (x — a)^ doubles the exponents at x = a, and likewise, a field 
extension of degree 2 can divide the exponents in half. The exponents at x = 1 of 
L4_3 must be divided in half to match (modulo the integers) the exponents of L4 at 
X = 1/16 and x = 00. So this field extension must ramify at x = 1/16 and x = 00, 
and this tells us that the field extension must be C(x) C C(x, \/l — 16x). We can 
write this latter field as C(^) where ^ = ^/l — 16x. With a pullback in C(^), the 
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a; = 1 singularity of L4 3 can be matched with the x — 1/16 and a; = cx) singularities 
of L4. 

A necessary (but not sufficient) condition for a homomorphism between two 
operatorqa to exist, is that the exponents of the singularities must match modulo 
the integers, and the number of logarithms must match at every singularity. But once 
one knows that the puUback for L4 must be in C(^), it suddenly becomes easy to find 
a puUback that meets this necessary condition. Once the puUback is found, we can 
check if a homomorphism exists (and if so, find it) with DEtools[Homomorphisms] in 
Maple. 

8.2. The 4F3 result 

Seeking for 4i^3 hypergeometric functions up to homomorphisms, and assuming an 
algebraic pull-back with the square root extension, (1 — 16 • w^)^' ^, we actually founqjj 
that the solution of L4 can be expressed in terms of a selected 4i*3 which is precisely 
the Hadamard product of two elliptic functions 

.F,{[l,y-,l],[l,l,l];^) (112) 

- ^^i([i^Mi];^)-2F,([i,i],[i] 

where the pull-back z is nothing but s^ with x = w'^ , where w is the natural variable 
for the x'"^'s n-fold integrals 8, 33l, w = s/{2 (1 -I- s^)): 

Let us recall that t = k^ = s'' and that EllipticK{k) in Maple is the integral 
with a fc^ in the square root, so t^ ~ k'^ — s^ 

^■EllipttcKiy) = 2i^i([i,i],[l]; y2^ (114) 

therefore the solution of L^ is expressed in terms of the Hadamard square of 
EllipticK , yielding in Maple notations: 
/I 1 1 1 
4^Hl2'2'2'2]'f^'^'^]' 

= 2Fi([i, i],[l]; t^)^^F,([\,\],[l]; e' 



2 2 

= - • EllipticK{t) • - • EllipticK{t), (115) 

extremely similar to the previously seen Hadamard product (|100p . 

8.3. The Calabi-Yau result 

This result could be seen as already achieving the connection with elliptic curves we 
were seeking for. In fact, looking at the Calabi-Yau list of fourth order operators 
obtained by Almkvist et al. [TT], one discovers that this selected 4i^3 hypergeometric 
function actually corresponds to a Calabi- Yau ODE. This is Calabi- Yau ODE number 

§ Here the two operators are: L4 and a puUback of -L4,3, both viewed as elements of C(5)[Dj] 
f Details will be given in forthcoming publications. 
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3 in page 10 the Almkvist et al. list, (see Table A of Calabi- Yau equations page 10 
in[n]). 

Remark 1: The hypergeometric function (|112p also corresponds to the hyper 
body-centered cubic lattice Green function [73 ( 178 ) : 

F(0,z) = (116) 

dkidk2dk^dk4^ 



/o Jo Jo Jo 1 ~ z ■ cos{ki) 003(^2) cos(/c3) cos(fc4)' 

It may well be, following the ideas of Christoll \79\ I^OJ /. that the Calabi- Yau three- 
fold corresponding to (I116p . (11121) similar to (J104I) . is nothing but the denominator of 
the integrand of (|116l) . (1 — z-cos(fci) cos(fc2) 003(^3) 003(^4)), written in an algebraic 
way {zi = exp(ifci)): 

^f4\\\,\ui.i.i\,^) - (117) 

dzidz2dz^dz4^ 



! zi Z2 z-i Z4 - (1 + z{) ■ (1 + zi) ■ (1 + zf) • (1 + z|) • 2 ■ 

Along this line 

8 zi Z2 zj, Zi - (1 + zl) ■ (1 + z|) • (1 + zl) .{I J^ zl)-z = 0, 

is a genus-one curve in (zi, Z2) which has to be seen on the same footing as the 
three- fold Calabi- Yau manifold (|104l) . 

Remark 2: Note that the series expansion of pi2l) for the inverse 1/ z of the 
pull-back, is a series with integer coefficients ir|j] the variable w: 

.^,([1 11 i[M.l'- ''-<'- '"■'"">"" 



2' 2' 2' 2^'^ ' ' ^' Vl + (1 -16-u;2)i/2. 
= 1 +16u;^ +512u;i° + 11264 w^^ ^ 212992 lu^^^ 
+ 3728656 w^^ -h 62473216 w^^ + 1019222016 v?^ (118) 

+ 16350019584^22 + 259416207616^2^ 
+ 4086140395520^26 + ... 

We also have this integrality property for 256 z : 

4-F3([^, ^, ^, ^], [1,1,1]; 256 z) = 1 +16 z +1296 z2 +160000^3 

+ 24010000 z^ + 4032758016 z^ + 728933458176 z^ 

+ 138735983333376 z'^ + 27435582641610000 z* 

+ 5588044012339360000 z^ + 1165183173971324375296 z^^ 

+ 247639903129149250277376 z" (119) 

+ 53472066459540320483696896 zi2 ^ .., 

The solution of L4, had been seen to be a series with integer coefficients (see (|92|) ). 
Now that we know that L4 has a Calabi- Yau interpretation, this integrality property 
can be seen as associated with mirror maps and mirror symmetries (see section (|9]) 
below) , as well as inherited from the Hadamard square structure (see ( [Appendix A.l[ ) 
below) . 

t This is not true in z or s. 



Ising model, modular form,s and Calabi- Yau 29 

8.4- Speculations: 4-F3 generalizations and beyond 

Let us consider a few generalizations of (|112p . the selected 4F3 we discovered for L4. 
More generally, the hypergeometric function 

^ /.I 1 1 1 , ^ 1 \ 

4^3 (^[2 + 9, 2 + ''' 2 ^ *' 2 ^ ^^'^"' ^ "^' "^ ^"^'^ ^ "^'' ^J' 
corresponds to the linear differential operator: 



7: + s] ■ (e + - + t 



^n. 


m, p: q, r, 


s.t — 


{' 


1 
+ 2 + 


,).(. 




1 

X 


(e + n 



2 / V 2 

(120) 

For any integer n, m, p, q, r, s, t all these operators (|120p are actually equivalent 
(homomorphic, see ( [Appendix B[ )). Therefore, all these linear differential operators 
p20p are homomorphic to (|120p for n = m = p = q = r = s = t = 0, which is 
actually a Calabi- Yau equation. 

We have also encountered another kind of generalization of (|112p : the Hadamard 
powers generalizations (jA.4p (see ( [Appendix A.l[ )). 

It is thus quite natural to consider the linear differential operators corresponding 
to generalizations like 

nFn-l {[^ + Pi, ■■■,^ + Pn], [1 + Ql, • • • , 1 + qn-l]] 4" • 2; j , 

where the p^'s and gi's are integers and see, if up to homomorphisms and rational 
or algebraic^ pull-backs one can try to understand the remaining quite large 
order operators L12 and L23, in such a large enough framework. There is no 
conceptual obstruction to such calculations. The obstruction is just the "size" of 
the corresponding massive computer calculations necessary to achieve this goal. 

9. Mirror maps for the Calabi- Yau 4i^3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; 256a;) 

An irreducible linear differential equation is said to be of Maximal Unipotent 
Monodromy (MUM) if all the exponents at are zero (one Jordan block). This is 
the case for all the hypergeometric functions 



^^«-l([^, ^, ■■•, ^],[i, 1, •••, i];4^ 



The hypergeometric function 

iF^{\^,\,\,\],[l,lA]; 256.t), 
which is MUM, corresponds to the fourth order linear operator 

x'^ ■{l-2hQx)-Dl + 2 2)3 -(3- 1024 x)-!)^ 

+ a;2.(7-3712a;)-L>2 + x • (1 - 1280a;) • D^^ -16a;, 
or, using the homogeneous derivative 9: 

6^ ~256-x-(e + -Y, e=x--^. (121) 

V 2/ dx 

^ Not too involved in a first approach, just square roots extensions. 



Ising model, modular form,s and Calabi- Yau 30 

It verifies the symplectic condition [S5]: 

1 1 3 da2 3 daa 1 d'^a^ 

a, = --a,-as - ^ ■ a, +— - ^ • «3 • ^ " 2 ' ^^' ^^^'^ 

for the monic order- four operator: D^ + a^ ■ D^ + 02 • D^ + ai ■ D^ + oq. 

Condition (|122p is nothing but the condition for the vanishing of the head 
coefficient of D^ of this exterior square (see Proposition 3 of 62 ). The exterior 
square of (|12ip is an irreducible order-five operator, instead of the order-six operator 
one expects at first sight. 

This opens room for a non- degenerate alternate 2- form invariant by the 
(symplectic) group 5P(4, C). Actually (|12ip has a 5P(4, C) differential Galois group. 

Remark: Let us compare this situation with the one for a "similar" order-four 
operator homomorphic to (jl21[) : 

9^ - 256 • a; • fe* - i V (123) 

The exterior square of (jl23p is an order-six operator which is the direct sum of an 
order-five operator homomorphic to the order-five exterior square of \121]) and an 
order-one operator 

^- +n oL , (124) 

(1 — 256 x) • X 

which has the simple rational solution (1 — 256a:)/a;. 

These two operators, (|12ip and (|123l) . have the exact same SP{4:, C) differential 

Galois group^ but, nevertheless, their corresponding exterior squares do not have the 

same orderly. The situation for (|12ip can be thought as an "evanescence" of the 

rational solution. 

9.1. Mirror maps in a MUM framework 
The solutions of (|12ip read: 



yo 



1 + 16 a; -I- 1296 x"^ + 160000 x^ + 24010000 x^ + 4032758016 x^ 



728933458176 a;^ + •••, 



Vi = Vo •ln(x) + yi, with: 

^,,„ 2 2368000 . 365638000 , 104147576064 , 

2/1 = QAx + 6048 x^ H x'^ ^ x^ H x^ 

3 3 5 

19045884743424 <, 25588111188676608 , 
+ 1 ^ + 35 ^ +••• 

h\{xf 
y2 = yo ^ h yi ■ \-a{x) + 2/2, with: 

, 8182400 . 1374099650 , 
y2 = 32 X + 5832 x^ -h x^ -f x^ 

685097536032 . 129379065232032 . 

H x^ H a;'" -f • • • 

25 25 

ln(a;)^ , ln(a;)^ _ , , , 

2/3 = yo 7 ^ 2/1 o ^ y2 • ln(a;) + yg, 

D 2 

_„^ n 10334080 , 1110845155 , 

W3 = - 64 a; - 4296 a;^ — x^ — x^ + ■■■ 

yi 27 27 

§ Since the log-degree of these operators is equal to four, the order of these exterior squares is at 
least five. 
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Introducing the nome (321 ^U\ q: 

q^ exp(—) = x-expf — V (125) 

one finds the expansion (with integer coefficients) of the nome q: 

q = X +Mx^ +7072a;^ + 991232 a;'' + 158784976 a;^ 

+ 27706373120 x^ + 64 x^ + 5130309889536 x^ 

+ 992321852604416 x^ 

+ 198452570147492456 x^ + 40747727123371117056 x^° 

+ 8546896113440681326848 x" + 1824550864289064534212608x^2 

+ 395291475348616441757137536a;" 

+ 86723581205125308226931367936 x^^ 

+ 19233461618939530038756686458880 x" 

+ 4305933457394032994320115176046592 x^^ 

+ 972002126960220578680860300103711764 x^^ 

+ 221026060926103071799983313019509871872 x" +••• (126) 

and, conversely, the mirror map |391 1401 162] reads the foUowing series with integer 
coefficients (x = z{q{x))): 

z{q) = q -64^2 +ll2Qq^ -38912^4 - 1536464 g^ 

- 177833984 g*^ - 19069001216 g^ - 2183489257472 g^ 

- 260277863245160?^ - 32040256686713856 g^" 

- 4047287910219320576 g" 

- 522186970689013088256 gi2 

- 68573970045596462152576 g" 

- 9140875458960295169327104 g" 

- 1234198194801672701733531648 q^^ 

- 168503147864931724540942221312 q^^ 

- 23230205873245591254063032928212 q^"^ 

- 3230146419442584387013916457526784 g" +••• (127) 
The Yukawa coupling [39|, 162) 



has the following (integer coefficients) series expansion: 
K{q) = 1 +32g + 4896 g^ + 702464 g^ 
+ 102820640 q^ + 15296748032 g^ 
+ 2302235670528 g^ + 349438855544832 g^ 
+ 53378019187206944 g^ + 8194222260681725696 g^ 
+ 1262906124008518928896 g" + 195269267971549608656896 g" 
+ 30273112887215918307768320 g^^ 
+ 4703886698057200436126953472 g" 



92= 77 -JT-r^ ocg ■,',2 ' (1^1) 
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+ 732300206865552210649383895040 g^^ 

+ 1 14192897568357606610746318782464 gi5 

+ 17832557144166657247747889907477280 gi^ (129) 

+ 2788280197510341680209147877101177216 g^^ 

+ 436459641692984506336508940737030913792 gi^ + ••• 

The nome series (|126p has a radius of convergence R — 1/256, corresponding 
to the z = 1/256 singularity. The m,irror map (jl27p . as well as the Yukawa series 
(|129p . have a radiufjjl of convergence R ~ 0.0062794754 • • • , corresponding to the 
singularity Qs — 0.0062794754 • • • given by: 

Qs = exp( — I, with: (130) 

a;o = 4F3[[-, -, -, -],[^A,'^]-A), and: 

OO 

xi = ^ 4 • r(n + 1/2)"^ • (*(n + 1/2) - ^{n + l))/r{n + 1)^/tt^. 

Introducing the rational functioqjl 

1 327680^2 _i792z + 5 

2 z2- (l-256z)2 

one finds that the mirror map (jl27p actually verifies a generalization of (|42p . namely 
the so-called quantum deformation of the Schwarzian equation (see (4.24) in [39j): 

f.=«M..i^ 14-i-(S^ 

where the derivatives are with respect to r, the ratio of the first two solutions 
T = j/i/j/o (log of the nome (|125p ). and {z, t} denotes the Schwarzian derivative 
63]). 

The rhs of (jl32p generalizes the very simple rational function rhs we had on the 
Schwarzian equation P^ (see the related footnote). The rhs of (I132p is basically 
a transcendental function depending on the Yukawa coupling function (jl28p . It is 
natural to try to obtain a (non-linear) ODE bearing only on the mirror map z{q), and 
not the Yukawa coupling function (|128p as well. This can be done (see ( [Appendix CP ) 
with a (complexity) price to pay, namely that these higher order Schwarzian (non- 
linear) ODEs are much more involved ODEs of much larger order. 

9.2. Higher order Schwarzian ODEs on the mirror map 

Actually, we have also obtained the higher order Schwarzian (non-linear) ODE (see 
(4.20) in [39j). verified by the mirror map (|127l) . It is an order-seven non-linear ODE 
given by the vanishing of a polynomial with integer coefficients in z, z', z", • • • , z"', 
having 1211 monomials in z, z', z", • • • ,z^^'. This polynomial of degree 18 in z, 24 
in z', twelve in z", six in z^^\ four in z^^\ three in z^^\ two in z'^' and one in 
z^^-' can be downloaded in [83 to check that (I127P actually verifies this higher order 
Schwarzian (non-linear) ODE. 

^ Along this line of radius of convergence of the mirror map and related Schneider-Lang transcendence 
criteria, see 181]. 

X Which corresponds to take the values ro = 16, r2 = 384 = (3/2) ■ 256, rs = 512, r4 = 256, 
in |39]. 



Ising model, modular form,s and Calabi- Yau 33 

One can verify that these higher order Schwarzian ODEs, on the mirror 
map, are actually compatible with the (renormahzation group, isogenics [33], ...) 
transformations q ^f q^ , for any integer n. Changing q — > q"^ in the mirror map 

z{q^) = g" -64g2" + 1120 g3" -38912 g4" _ 1536464 g^" 

- 177833984 g^" - 19069001216 9^" - 2183489257472 g^" 

- 260277863245160 g^" + ••• (133) 

one finds that this new function is still a solution of the higher order Schwarzian ODE. 
Conversely, one can consider the reciprocal higher order Schwarzian ODE bearing 
on the log of the nome, t = \n{q) = t(z), seen as a function of z. It is an order-seven 
non-linear ODE given by the sum of 602 monomial terms: 

0= z« -(1- 256 • z)2.p,(z, r',r", •••r(-^))-^^'^ + ■•• 
+ (l-256-z)-Q,(z, T,---,r(''),r(^)) 

+ i?,(z, r, T',---,T(5))-r', (134) 

where the t^^^'s are the m,-th z-derivative of t(z), and where Pr and Rr are 
polynomials of the z and the r'^™^'s derivatives (see (|D.3p below). 

Let us consider the Moebius transformation (homographic transformation) on the 
log of the nome 

T — > — ., (135 

c- T + a 

which transforms, as far as the ^-derivatives are concerned, in an increasingly involved 
way with increasing orders of derivation: 



ad — cb 



2 



[c-T +d) 

{ad — cb) „ (ad — cb) ■ c ,2 

2- -T^ ■ T , 



(136) 



(c • T -f d)2 (c-T + d)^ 

It is a straightforward calculation to verify that the higher order Schwarzian ODE 
(jl34p is actually invariant by the Moebius transformation (jl35p and its deduced 
transformations on the derivatives (J136I) . Do note that we do not impose ad ~cb = 1: 
we are in GL{2, Z) not in SL{2, Z). The previous symmetry q — > g" (see (|133p ) of 
the higher order Schwarzian ODE corresponded to r — >■ n-r. We have here a GL{2, Z) 
symmetry group of the higher order Schwarzian ODE, corresponding to the extension 
of the well-known modular group SL(2, Z) by the isogenics (exact representation of 
the renormahzation group 25!) r -^ n-r, which extend quite naturally the modular 
group and isogenics symmetries encountered with elliptic curves [25' . Even leaving the 
elliptic curves or modular forms framework, for some natural generalizations (Calabi- 
Yau are natural generalizations of Hauptmoduls) it was crucial to get mathematical 
structures with canonical exact representation of the renormahzation group 25 . 

10. Late comments: the integrability behind the mirror. 

It is beyond the scope of this very down-to-earth paper to give a mathematical 
definition of mirror symmetries, since mathematicians are still seeking for the proper 
general framework to define them (mixed Hodge structures, flat connection underlying 



Ising model, modular form,s and Calabi- Yau 34 

a variation of Hodge structuretp in the Calabi- Yau case [HHIHSIIHS]! T-duality [87 , toric 
frameworks like in Batyrev's construction of mirror symmetry between hypersurfaces 
of toric Fanq^ varieties [551, algebraic Gauss-Manin connections 'SDI, ...). 

More familiar to physicists, in particular specialists of integrability, is the notion of 
Picard-Fuchs [91] equatiorl^. Along this Picard-Fuchs line, the occurrence of Painleve 
VI equations is well-known for the Ising model [13l [T4| (see the Garnier or Schlesinger 
systems [93]). This Picard-Fuchs notion is central in any "intuitive" understanding of 
mirror maps [94] and other Calabi- Yau manifolds (namely compact Kahler manifolds 
with Ricci-flat Kahler metrics). 

As far as the "proper integrable framework" of this paper is concerned let us 
underline the following comments. The 2-D Ising model is a well-known free-fermion 
model with an elliptic parametrization. This elliptic parametrization is, of course, 
a straight consequence of the Yang-Baxter integrability of the model (here the star- 
triangle relation), and, therefore, it is not a surprise to see elliptic functions in the 
integrals of correlation functions (see for instance [531111117]). Along this line, even 
Painleve VI equations can be seen as a Gauss-Manin deformation of an elliptic function 
second order ODE [Ml dH [9H] . 

However, it is crucial to note that the elliptic parametrization is not one-to- 
one related with a Yang-Baxter integrability: the sixteen vertex model, which, 
in general, is not Yang-Baxter integrable, has a canonical (compulsary !) elliptic 
parametrization [99], the elliptic parametrization being a consequence of the 
integrability of the birational symmetries |23j of the model, and of course not of a 
Yang-Baxter integrability that does not exist generically for that very model. The 
free-fermion character of the square Ising model is of course crucial in Wu et al. [2] 
(Pfafhan, Toeplitz, ...) calculations to write explicitly the ^("^'s as integrals of some 
integrand algebraic in some well-suited variables. However, Guttmann's paper |74| 
makes crystal clear with miscellaneous examples of Green functions for many lattice 
statistical mechanics, or enumerative combinatorics, problems in arbitrary lattice 
dimensions, that Calabi- Yau ODEs emerge in a lattice statistical mechanics framework 
which is (at first sight) quite remote from Yang-Baxter (tetrahedron, ...) integrability, 
and even more from free-fermion integrability. 

If the occurrence of linear differential operators associated with elliptic curves for 
square Ising correlation functions, or form factors, is not a surprise [77], the kind 
of linear differential operators that should emerge in quite involved highly composite 
objects like the n-particle components x^"-* of the susceptibility of the square Ising 
model, is far from clear. We just had a prejudice that they should be "special" 
and could possibly be associated with elliptic curves. In fact, even if Yang-Baxter 
structures were the "deus ex machina" behind these "special" linear differential 
operators it seems impossible to use that property in any explicit calculation. 

Moving away from Yang-Baxter integrability to other concepts of "integrability" , 
we are actually using the following ingredients: we have n-fold integrals of an 
integrand which is algebraic in the variables of integration and in the other remaining 
variables. This algebraicity is the crucial point. As a consequence we know that 

X See 1821 for the introduction of the notion of variation of Hodge structures with their Gauss-Manin 

connections. 

§ Mirror symmetry, in a class of models of toric varieties with zero first Chern class Calabi- Yau 

manifolds and positive first Chern class (Fano varieties) was proven by K. Hori and C. Vafa 1881 . 

f It is known that if a linear differential equation with coefficients in Q is of Picard-Fuchs type, then 

it also describes an abstract variation of Q-Hodge structures, and it is globally nilpotent [92]. 
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these n-fold integrals can be interpreted as "Periods" of algebraic varieties and verify 
globally nilpotent [5] linear differential equations: they are [331I351I25] "Derived From 
Geometr'i/\ However, inside this "Geometry" framework [lOOj (in the sense of the 
mathematicians) theoretical physicists are exploring^] "Special Geometries". These 
linear differential operators factorize in irreducible operators that are also necessarily 
globally nilpotent [9] . When one considers all the irreducible globally nilpotent linear 
differential operators of order TV, that we have encountered, (or the one's displayed 
by other authors in an enumerative combinatorics framework [73[ I74j , or in a more 
obvious Calabi-Yau framework [511 IZJ), it apperas that their differential Galois group 
are not the SL{N, C) or extensions of SL{N, C), groups one could expect generically, 
but selected SO{N), SPiN, C), G2, ... differential Galois group [TOT]. These are, 
typically, classification problems in algebraic geometry and/oi|fl differential geometry. 
Our linear differential operators are, in fact, "special" globally nilpotent operators (G- 
operators). This paper can be seen as an attempt, through a fundamental model, the 
Ising model, to try to characterise the "additional structures and properties" of these 
globally nilpotent operators. In the simple example of pFq generalised hypergeometric 
functions, only operators with nFn-i solutions [60] can be globally nilpotent. In a 
hypergeometric framework we are thus trying to see the emergence of "special" nFn-i 
hypergeometric functions. 

The last results, displayed in this paper, show clearly, with the SP{4, C) 
differential Galois group of L4 for x^^\ that these "special geometries" already 
emerge on the Ising model which, therefore, does not restrict to the theory of elliptic 
curves jl02j (and their associated elliptic functions and modular forms). To define 
these "special geometries" is still a work in progresijj but it seems to be close to 
concepts like the concept of modularity (for instance, integrality of series) and other 
mirror symmetries |103j . 

11. Conclusion 

All the massive calculations we have developped during several years [TJ [S] [6j [T] IH [21] 
on the square lattice Ising model give coherent exact results that do show a lot 
of remarkable and deep (algebraico-differential) structures. These structures all 
underline the deep connection between the analysis of the Ising model and the theory 
of elliptic functions (modular forms, selected hypergeometric functions |25j, modular 
curves, ...). In particular we have actually been able to understand almost all the 
factors obtained in the analysis of the x Sj ^^ linear differential operators ''associated 
with elliptic curves" . Some linear differential operators have a very straightforward 
relation with elliptic curves: they are homomorphic to symmetric powers of Le or 
Lpc, the second order operators corresponding to complete elliptic integrals E or K. 
We showed, in this paper, that the solutions of the second and third order operators 
Z2, F2, F3, L3 operators, can actually^ be interpreted as modular forms of the elliptic 

§ Sometimes, without knowing it, like Monsieur Jourdain (Le bougeois gentilhomme) . 

tl Given a classification problem in algebraic geometry, using the mirror duality one can translate 

it into a problem in differential equations, solve this problem and translate the result back into 

geometry. 

f For the next X ^''^i ™ ^ 7, we just know that the corresponding differential operator are globally 

nilpotent. We can only conjecture the emergence of these "special geometries", as we already 

conjectured the integrality of the ^'■"''s series in well-suited variables (see equation (8) in [5]). 

^ At first sight it looks like a simple problem that could be solved using utilities like the "kovacicsols" 

command |104| in maplelS: it is not even simple on a second order operator. 
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curve of the Ising model. These resuhs are aheady quite a "tour-de-force" and their 
generalization to the much larger (and involved) operators L12 and L23, seems out 
of reach for some time. The understanding of the "very nature" of the globally 
nilpotent fourth order operator L4 was, thus, clearly a very important challenge 
to really understand the mathematical nature of the Ising model. This has been 
achieved with the emergence of a Calabi- Yau equation, corresponding to a selected 
4i^3 hypergeometric function, which can also be seen as a Hadamard product of the 
complete elliptic integrals K, with a remarkably simple algebraic pull-back (square 
root extension (J113I) ). the corresponding Calabi- Yau fourth-order operator having a 
symplectic differential Galois group 5*^(4, C). The associated mirror maps and higher 
order Schwarzian ODEs present all the nice physical and mathematical ingredients 
we had with elliptic curves and modular forms, in particular an exact (isogenics) 
representation of the generators of the renormalization group, extending the modular 
group SL{2, Z) to a GL{2, Z) symmetry group. 

We are extremely close to achieve our journey "from Onsager to Wiles" (and now 
Calabi- Yau ...), where we will, finally, be able to say that the Ising model is nothing but 
the theory of elliptic curves, modular forms and other mirror maps and Calabi- Yau. 
Do note that all the ideas, displayed here, are not specific of the Ising mode^ and can 
be generalized to most of the problems occurring in exact lattice statistical mechanics, 
enumerative combinatorics [731 IZE]) particle physics, ..., (the elliptic curves being 
replaced by more general algebraic varieties, and the Hauptmoduls being replaced 
by the corresponding mirror symmetries generalizations |50j). 

We do hope that these ideas will, eventually yield the emergence of a new 
Algebraic Statistical Mechanics classifying all the problems of theoretical phvsic affl 
on a completely (effective) algebraic geometry basis. 
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Ising model, modular form,s and Calabi- Yau 37 

Appendix A. Hadamard product of operators depends on the expansion 
point 

The Hadamard product of two operators depends on the point around which the series 
expansions are performed, and, hence, the Hadamard product of the series-solutions. 
Let us consider the complete elliptic integral K which is already a Hadamard 
product 

-■EllipticK{4:-x^/'^) = (1 -4a;)"i/2^(l -4a;)"i/^ (A.l) 

TT 

and the order-two linear differential operator for EllipticK{x^/^) 

m^'l^-n.-l—L^. (A.2) 

[I ~ x) ■ X 4: [1 ~ X) ■ X 

The Hadamard square of a linear differential operator does depend on the 
point around which the series are performed. For EllipticK with its three 
singularities 0, 1, cxd, the Hadamard square yields the same operator of order- four 
for the three expansion points 0, 1, oo, this order four operator corresponding to 

4^3 ([5, 5, 5, 5]Jl,l,l]; 2 j, namely: 

, (3 -4x) . 1 (14 -29 a;) ^ 

Dl -K 2 j- -^ -01+- )- ^ ^ • Dl 

^ {l-x)-x ^ 2 (1 - x) • x2 -^ 

^ {l-x)-x^ " 16 (1 - a;) • x3 ' ^ ^ ^ 

for the three expansion points 0, 1, 00. However, for a generic expansion point, x = c 
(where c 7^ 0, 1/2, 1, cx)) one gets an order-six linear differential operator. 

Appendix A.l. Hadamard powers generalizations 

Let us denote HadJ^lF) — F -k F -k ■ ■ ■ -k F, the Hadamard product of F, n-th time 
with itself. Relation (|112l) can straightforwardly be generalized to arbitrary Hadamard 
powers of the complete elliptic integral K, which, as we know, plays a crucial role in 
our analysis of the Ising model [T31 E] : 

^, •••,^],[i, •••,1]; le"-^) (A.4) 

4 1, 

These hypergeometric functions |108j are solutions of the 2n-th order linear differential 
operator: 

6)2" _ 16". z. [0 + -j . (A.5) 

These relations are a subcase of the (slightly) more general Hadamard power relations 

1 



i?arf"(/i([^, ^],[1]; 16- 



,F„_i([i, ...,i],[l, ...,1]; 4'^ 



Had'- 



these last hypergeometric functions being solutions of the n-th order linear differential 
operator: 



4" 



(. + -). (A6) 
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Do note that the corresponding („)F(„_i) hypergeometric series are, actually, series 
with integer coefficients: 

1 + 2" • z + 6" • z^ + 20" • z^ + 70" • z'^ + 252" • z^ + 924" • z^ 
+ 3432"-z^ + 12870" -z*^ + 48620" • z^ + ••• (A.7) 



Appendix B. Equivalence of the Jfc„; Equivalence of the Q, 



n, ra, p\ q,r, s,t 

• All the operators Jk,n defined by (|106p are homomorphic. This can be seen 
recursively from the two operator equivalences 

Ul ■ Jk,n — Jk,n+1 ' W, Vl • Jk,n = Jk+l,n ' V, 

where 

4 n 
U = (1-16 a;)- 6*3 

X 

(16 (8 n + 6 /en + /c^ + 5 n^) • a; - ik + ^n) (k + n)) ■ 9"^ 

X 

- (16 (fc + n + 1) (/c^ + 3 fcn + 5 n + 2 7i2) • 6* (B.l) 

k 
+ - • {bn^ + 2kn + k^)) • 6* - 4 • Ik^ + 2kn + An + n^) [k + n + if , 

X 

4 n 
Ui ^ U + (1 -32 a;)- 6*2 

X 

(8 (4 fcn + 8 n + 3 n^ + /^2^ ^^^yg 

X 

- 16- (fc + n + 2)(fc2 + 2fcn + 4n + n2) (B.2) 

and where V (resp. Vi) is lA (resp. lAi) where k and n have been permuted. 

• Let us now show that the order- four linear differential operators ^n.m,p;q,r,s.tj 
corresponding to the hypergeometric functions (jl20[) . are homomorphic. 
One has: 

^ ^n, m, p;q,r, s,i'lt'i^^iJ^I \y^\'n-\-L\-\ ln-\-l. ni. p: q. r.s.t^ 

^n,m,p;q, r, s,t'(^ + "*+ll = l9 + ni+l\- iln, m+1, p; q, r, s,t: 
^'n,m,p;g,r, s,i'lt'iJ'i-L| [^iPi^j'^^ji^m, p+l; q,r, s, t- 

and: 

Q O/ n^ m^ p; q^ r^ s^ t ^^n,?n,p; g+1, r, s, t ' [^ ' Q ' "^ 

\u -\- r 17^1 ' ^hi.ra^p] q^r^ s^t — ^'n,?7i,p;g, r+1, s,i'lt' + ^i~), 
I" ^" ^ ^j ^^' ^^' P! 'j'l ^1 ^1 ^ ^^^ ^*^' P; ^5 ^5 s+1, i' If^ ^" ^ ^/' 

A simple composition of all these elementary relations show, by recursion, that the 
^n,m,p; 9, r, s, t's are all homomorphic. 
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Appendix C. Getting the higher order Schwarzian ODEs 

Let us call £4 the order four linear differential operator corresponding to the Fuchsian 
ODE of the hypergeometric function 

4F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; 256a;). (C.l) 

The formal solutions of £4 are denoted j/3, y2, j/i and yo, where the subscript is for 
the higher exponent of the log. The nome map reads (see (I126p ) 

qix) = expf— ) = X +64a;2 + 7072 a;^ + 9912320;^ + ••• 

Our aim is to obtain the non linear ODE of q{x). 

Considering the log nome map as yi/yo = ln((?), and, differentiating both sides, 
gives 

y'l yo -yiy'o = yl- — (C.2) 

where the left-hand side is the wronskian of the solutions j/i and yo- This wronskian 
is solution of an order- five linear differential operator £5 , which is the exterior square 
of £4. We have then three equations 

A(yo--) = 0, £4(2/0 -Mg)) = 0, £4(2/0)= 0, (C.3) 

to which we add the following two equations obtained by differentiating the last two 
equations 

D, ■ £4 (2/0 • Hq)) =0, D,- £4 (2/0) = 0. (C.4) 

We solve this system of five equations in the unknowns i^'"^ — d'^yo/dx'^, to obtain 

F(") = B„-^(°), n= 1, 2, •••,5 (C.5) 

The i3n's depend on x, q{x) and its derivatives up to seven. All these relations should 
be compatible, by which it is meant that the derivative of one relation with n gives 
the relation with n + 1. The result is then 

— Bn +Bn-Bi -Bn+1 =0, 71= 1,2, •••,4. (C.6) 

The gcd of these four nonlinear differential equation is the ODE oiq{x). The nonlinear 
ODE of q{x) contains the derivatives q'{x), ■ ■ ■ , (/'■^^(x) with degrees, respectively, 12, 
16, 8, 5, 4, 3, 2 and 1. The nonlinear differential equation of the mirror map (see 

dMD) 

X{q) ^ q - 64 g2 +1120 9'"^ -38912 g^ -1536464 g^ +••• 

can be obtained by using X {q{x)) = x and the ODE of q{x). The nonlinear ODE of 
X{q) involves the derivatives A(°'(g), ^(^^(q), • • • ,X^'^^q) with degrees, respectively, 
18, 24, 12, 6, 4, 3, 2 and 1. 

Appendix D. Higher order Schwarzian ODEs for the hypergeometric 
function 4i^3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; 256 z) 

Actually the non-linear ODE on the mirror map (|127p is of the form 

= (1 - 256 • z) • S{z, z', z" ■■■ z(^), z^^\ z^^)) 
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+ z'i2.p(z, z',--- z(5))^ (D_l) 

F(z', z ••• 2(5)) = -18z'3.(188416z'4 +6z'z(3) -9z"2).z(5) 
+ 1909301941633024 z'i2 _ 695516266496 z'^z^^^ 
+ 12582912 ^82912 z"^ - 345 z^^A ■ z'^ - 4 (z^^))^) z'^ 

+ 26046627840 z'^ z" z^^) + 2430 z' z"^ z^^) 

- 33914880 (768 z"3 - (z^^))^ - z^^^ z") • z'^ 

- 203489280 z'^z"2z(3) + 135 (1130496 z"* + (z^''))^) z'* 

- 180 z(3) . (3 z(4) z" - 1620 z'2 z"2 (z^^))^ _ 1215 z"^. 

This is a non-linear ODE where the derivatives z*^"-* are the n-th derivatives in 
T = ln(g) of the mirror map z{q). 

Conversely, the ODE on the function r = ln(g) = t{z) is an order-seven non- 
linear ODE given by the sum of 602 monomial terms: 

0= z*^ -(1- 256 • z)2. P,(z, T',r", •••r(5))-T(^) -)-••• 

+ (1 - 256 •z)-Q,(z,r,---,T(«\T(")) 

+ i?.(z, r, r',...,T(5)).r', (D.2) 



where: 



i?,(T, t',---,t(-^)) = 18 • (9 t"2 -6 t(3)t' + 188416 r'2).T(5) 
- 1043274399744 t"V -I- 16957440 • (256 r' - t") • r' • r^^^ 
-h 13023313920 r"3 -I- 1909301941633024 t'^ -I- 360 (r^^')^ 

+ 135 T ■ (t^*) f + 4 • (^173879066624 t'^ - 135 r" t'*) 
- 4341 104640 t" r' + 4239360 r"^^ • t^^^ . 



The non-linear ODE (ID.2P is a homogeneous polynomial expression of degree four in 
the seven derivatives (r', r", • • • t^'''). 

Rewritten in g(z) the non-linear ODE is an order seven non-linear ODE given by 
the sum of 2471 monomial terms in (z, g, q' q" ■ ■ ■ g^^-*): 

= z^ • g^ • (1 - 256 • zf ■ Pg{z, q\ q" , ■ ■ ■ g^^^) • g^^^ -^ • • • 

+ (l-256-z).Q,(z, 9,...,g(6)^^(7)) (D,3) 

+ R,iz,q,q',---,q^'^)-q'-q\ 

where Pq and Qq are polynomials of z, q', q", ■ ■ ■ g^^^ and z, q' , q", ■ ■ ■ g^^^ 
respectively, and where Rq : 

i?,(z, g, g',...,g(5)) = 

18g4- (l88416gV' -6g'g(3)^2 +9 ^"2^2 _3^'4^ .^(5) 

+ 4341104640- (g'2g6g(4) + g*^" g'* - q^q'^) 

- 1043274399744 g'g*^g"2 - 270g'^g(*) q^ + 135 [q'-^'^fq^q' 

+ 695516266496 g'V?^^^ + 13023313920 gV'^ 

-I- 347758133248 g'V + 1909301941633024 g'V 

-I- 13565952 g'V - 17364418560 g'g^g^^^ g" - 36 g'^ 

+ 1215 g"^ g' g^ + 648 g"^ g'^ g^ - 432 g(^) g'^^g^ + 360 g^^'^ g*^ 
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+ W20{q^^'>q"q\"' - q"^q'^q^ ~ q^^W'^q^q'^) 
+ 540 -(g'^ 9" 9^9(4) - g(3)g"^6^(4)) 

+ 16957440 • (q^ g"2 g^^) _ ^' ^6 ^(4) ^/, 

One can verify that q — Constant is a solution of (jD.Sp . Furthermore the 
expansion (jl26p multipHed by an arbitrary constant is stih a solution of (|D.3|) : 

q ^ Co-{z + 6Az^ + 7072 z^ + 991232 z^ + 158784976 z^ + • • •) 

which is natural since (|D.3p is a linear ODE on derivatives of ln(5). 

Remark. Such non-linear ODE is a "machine" to build series with integer 
coefhcients. For instance, if we explore the solutions of (ID.3|) . of the form q = 
z^ + ■ ■ ■ , one gets: 

q = z^ + 128 z^ +18240 z^ +2887680 z^ +494460832 z^ 

+ 89757208576 z^+ 17035431116800 z^ + 3347987811139584 z^ 

+ 676624996390235600 z^° + 139902149755519715328 z" 

+ 29480532176870291252224 z^2 + 6312281522697932105646080 z" 

+ 1370123593804106822389706240 z" 

+ 300913725420989219840662110208 z^^ 

+ 66768780541145654061810373951488 z^^ + ••• (D.4) 

In fact, this new series is nothing but the square of ()126p . Similarly one easily verifies 
that the cube of p^ 

q = z^ +192z^ + 33504 z^ + 5951488 z^ 

+ 1093928304 z'^ + 207935296512 z* 

+ 40712043902464 z^ + 8176029744758784 z^" + • • • (D.5) 

is also a solution of (jD.3|) , and this is also true for negative powers of the expansion 
(|126p . for instance the inverse (in the sense of the multiplication), 

q= - -64 -2976 z -348160 z^- 52017616 z^ -8802913280 z" 

z 

- 1608195557888 z^ - 309505032060928 z^ + • • • (D.6) 
is also a solution of (|D.3|) . and similarly, (jl26p to the power (—2) 

q= J- +— -1856- 315392 z- 50614176 z^- 8875323392 z^ 

z z 

- 1658793979904 z'' + ••• (D.7) 

is also a solution of (|D.3[) . 

This remarkable property, expression of the renormalization growpj^, is in fact 
a straight consequence of the fact that the non-linear ODE (jD.21) is a homogeneous 
polynomial of degree four in the seven derivatives (t' , t" , • • • r^'^' ) . 

Conversely, in the "mirror" , in the non-linear ODE (|D.ip , one can change q into 
A ■ q since the derivatives are all log-derivatives of q 

z = A-q -M- A^ -q^ + 1120 ■ A^ ■ q^ - 38912 ■ A'^ ■ q^ 

- 1536464 •^^- q^ + ••• (D.8) 

I We have an exact representation |25| of t — > Af ■ t or q — f q^ . 
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is also solution of (|D.ip . and, of course, one can easily check the symmetry 
corresponding to change q — ^ g", the new mirror map series (jl33p 

z(q") = g2n „64g''" +1120<7^" -38912g^" +••• (D.9) 

being also solution of (jD.ip . 



References 

[1] S. Boukraa, A.J. Guttmann, S. Hassani, I. Jensen, J.-M. Maillard, B. Nickel and N. Zenine, 

Experimental mathematics on the magnetic susceptibility of the square lattice Ising model, J. 

Phys. A: Math. Theor. 41 (2008) 455202 (51 pp) and arXiv:0808.0763 
[2] T. T. Wu, B. M. McCoy, C. A. Tracy and E. Barouch, Spin-spin correlation functions for the 

two-dimensional Ising model: Exact theory in the scaling region, Phys. Rev. B 13 (1976) 

316-374 
[3] B. Nickel, On the singularity structure of the 2D Ising model susceptibility, J. Phys. A 32 (1999), 

no. 21, 3889-3906. 
[4] B. Nickel, Addendum to "On the singularity structure of the 2D Ising model susceptibility" , J. 

Phys. A 33 (2000), no. 8, 1693-1711 
[5] A. Bostan, S. Boukraa, A.J. Guttmann, S. Hassani, I. Jensen, J.-M. Maillard, and N. Zenine, 

High order Fuchsian equations for the square lattice Ising model: x , J- Phys. A 42: Math. 

Theor. (2009) 275209 (32 pp) and arXiv:0904.1601vl [math-ph] 
[6] B. Nickel, I. Jensen, S. Boukraa, A.J. Guttmann, S. Hassani, J.-M. Maillard, and N. Zenine, 

Square lattice Ising model x ODE in exact arithmetics, J. Phys. A 43: Math. Theor. 

(2010) 195205 (24 pp), and arXiv:1002.0161.'2 [math-ph] 
[7] S. Boukraa, S. Hassani, I. Jensen, J.-M. Maillard, and N. Zenine, High order Fuchsian equations 

for the square lattice Ising model: x'-^K J- Phys. A: Math. Theor. 43 (2010) 115201 (22 pp) 

and arXiv:0912.4968vl [math-ph] 
[8] N. Zenine, S. Boukraa, S. Hassani and J.-M. Maillard, The Fuchsian differential equation of 

the squ are Ising model x'^-* susceptibility, J. Phys. A: Math. Gen. 37 (2004) 9651-9668 and 

[ar"Xiv:m ath-ph/0407060 
[9] A. Bostan, S. Boukraa, S. Hassani, J.-M. Maillard, J. -A. Weil and N. Zenine, Globally nilpotent 

differential operators and the square Ising model, J. Phys. A: Math. Theor. 42 (2009) 125206 

(50pp) and'arXiv:0812.4931 
[10] M. A. Barkatou and E. Pfliigcl, On the equivalence problem of linear differential 

systems and its application for factoring completely reducible systems International 

Conference on Symbolic and Algebraic Computation, (1998), Proceedings of the 1998 

international symposium on Symbolic and algebraic computation, (ISSAC '98) pp. 268-275 

http: //portal. acm.org/citation.cfm?id=281633 
[11] Homomorphisms command in Maple's DEtools package, source code available at 

http: //www. math. fsu.edu/f^hoeij/files/Hom/ 
[12] M. F. Singer, Testing reducibility of Linear Differential Operators: a Group Theoretic 

perspective. Applicable Algebra in Engineering, Communication and Computing 7, (1996), 

77-104 
[13] S. Boukraa, S. Hassani, J-M. Maillard, B. M. McCoy, J-A. Weil, N. Zenine, Painleve versus 

Fuchs, J. Phys. A 39 (2006) 12245-12263 Special issue for the centenary of the publication 

of the Painleve VI equation http://arxiv.org/pdf/math-ph/0602010 v3 
[14] S. Boukraa, S. Hassani, J-M. Maillard, B. M. McCoy, J-A. Weil, N. Zenine, Fuchs versus 

Painleve, J. Phys. A 40 (2007), 12589-605 Special issue for the SIDE VII conference 

(Melbourne) |http://arxiv\or g/pdf/math-ph/0701014 v2 
[15] D.H. Bailey, D. Borwein, J.M. Borwein, R.E. Crandall: Hypergeometric Forms for Ising-class 

Integrals Experiment. Math. 16, Issue 3 (2007), 257-276. 
[16] R. S. Maier, On rationally parametrized modular equations, J. Ramanujan Math. Soc. 24 (2009); 

no. 1, 1-73, and (2006) arXiv:NT/0611041 v4. 
[17] A.J. Guttmann and T. Prellberg, Staircase polygons, elliptic integrals, Heun functions, and 

lattice Green functions, Phys. Rev. E 47, Rapid Comm. (1993), R2233-R2236. 
[18] R. Maier, Algebraic hypergeometric transformations of modular origin. Transactions of the 

American Mathematical Society, 359, 3859-3885, (2007). 
[19] R. Maier, On reducing the Heun equation to the hypergeometric equation, J. Differential 

Equations 213, (2005), no. 1, 171-203. 



Ising model, modular form,s and Calabi- Yau 43 

[20] J-M. Maillard and S. Boukraa, Modular invariance in lattice statistical mechanics, Annales de 

la Fondation Louis de Broglie, Volume 26, nuinero special, pp. 287-328 (2001) 
[21] S. Boukraa, S. Hassani, J.-M. Maillard and N. Zenine, Singularities of n-fold integrals of the 

Ising class and the theory of elliptic curves, J. Pliys. A: Math. Tlieor. 40 (2007) 11713-11748 

and arXiv.org/pdf/0706.3367 
[22] J-M. Maillard, Automorphisms of algebraic varieties and Yang-Baxter equations, Journ. Math. 

Phys. 27, (1986), 2776-2781 
[23] S. Boukraa and J-M. Maillard, Let's Baxterise, J. Stat. Phys. 102, (2001), 641-700, and : arXiv: 

hcp-th/0003212 
[24] M. Barcau, Isogeny covariant differential modular forms and the space of elliptic curves up to 

isogeny, Compos. Math. 137 (2003), 237-273. 
[25] A. Bostan, S. Boukraa, S. Hassani, J-M. Maillard, J-A. Weil, N. Zenine and N. Abarenkova, 

Renormalization, isogenics and rational symmetries of differential equations, Advances in 

Mathematical Physics Volume 2010 (2010), Article ID 941560, (44 pp) 

[^Xiv:0911.5466v2 [math-ph] 
[26] B. C. Berndt and H. H. Chan, Ramanujan and the Modular j -Invariant, Canad. Math. Bull. 42 

(1999) 427-440 
[27] R. B. Griffiths and P. A. Pearce, Position-Space Renormalization- Group Transformations: Some 

Proofs and Some Problems, Phys. Rev. Lett. 41, (1978) pp. 917-920. 
[28] R. B. Griffiths and P. A. Pearce, Mathematical Properties of Position-Space Renormalization- 

Group Transformations, J. Stat. Phys. 20, (1979) pp. 499-545. 
[29] A. C. D. van Enter, The Renormalization- Group peculiarities of Griffiths and Pearce: What 

have we learned ?, Proceedings 1998 Marseille lUPAP-Statphys satellite Colloquium on 

Mathematical Results in Statistical Mechanics, eds. S. Miracle-Sol, J. Ruiz and V. Zagrebnov, 

pp. 509-526. World Scientific, 1999 and arXiv: |cond-m ath/9810405vl (1998) 
[30] H. H. Chan and M.-L. Lang, Ramanujan's modular equations and Atkin-Lehner involutions, 

Israel Journal of Mathematics, 103, (1998) 1-16. 
[31] H. Cohn, Fricke's Two-Valued Modular Equations, Mathematics of Computation, 51, Number 

184, (1988), 787-807 
[32] S. Boukraa, S. Hassani, J.-M. Maillard and N. Zenine, Landau singularities and singularities 

of holonomic integrals of the Ising class, J. Phys. A: Math. Theor. 40 (2007) 2583—2614 and 

[ajXi v:math-ph/0701016 v2 
[33] N. Zenine, S. Boukraa, S. Hassani and J.-M. Maillard, Square lattice Ising model susceptibility: 

series expansion method, and differential equation for x"^, J- Phys. A: Math. Gen. 38 (2005) 

1875-1899 and arXiv:hep-th/0411051 
[34] N. Zenine, S. Boukraa, S. Hassani and J.-M. Maillard, Square lattice Ising model susceptibility: 

connection matrices and singular behaviour of x^^' o-i^d x^\ J- Phys. A: Math. Gen. 38 

(2005) 9439-9474 and arXiv:math-ph/0506065 
[35] P. Candelas, X. de la Ossa, P. Green and L. Parkes, A pair of Calabi-Yau manifolds as an 

exactly soluble superconformal theory, Nucl. Phys. B359, (1991), 21—74. 
[36] C. F. Doran, Picard-Fuchs Uniformization and Modularity of the Mirror Map, Comm. Math. 

Phys. 212, 625-647, (2000). 
[37] C. F. Doran, Picard-Fuchs Uniformization: Modularity of the Mirror Map and Mirror- 
Moonshine, CRM Proc. Lecture Notes, 24, Amer. Math. Soc. 257-281, Providence and 

larXiv:math/9812162vl, (2000). 
[38] C. Krattenthaler and T. Rivoal, On the Integrality of the Taylor Goefficients of Mirror Maps, 

Duke Math. Journ. 151, Number 2, (2010) 175-218, and ar Xiv: 0907. 2577^ 2 [math.N.T] and 

|irXiv:0907.2578v2 [math.N.T] 
[39] B.H. Lian and S.T. Yau, Arithmetic properties of mirror map and quantum coupling, Comm. 

Math.Phys. 176 pp.163-191 (1996) and arXiv/hepth/9411234v3 
[40] B.H. Lian and S-T. Yau, Mirror Maps, Modular Relations and Hypergeometric Series II, Nuclear 

Phys. B 46, Proceedings Suppl. Issues 1-3, (1996) 248-262 and arXiv: hepth/950753vl (1995) 
[41] I. Jensen, Some exact results for factors occuring in the order 29 differential operator L29 

annihilating the series 0(^) = x'^'' ~ X^^^ /'^ + x'^Vl20, 

^tp://www. ms.unimelb.edu.au/~iwan/ising/Chi5.htmll 



and http://www.ms.unimelb.edu.au/~iwan/ising/series/Chi5/Lll.txt 
[42] I. Jensen, Series and other data related to the six-particle contribution x*- ' to the square lattice 

Ising model susceptibility, (2009) http://www.ms.unimelb.edu.au/~iwan/ising/Chi6.html 

and http://www.ms.unimelb.edu.au/~iwan/ising/series/Chi6/L17.txt 
[43] J-A. Weil, 2^1 and 4F3 type solutions of some order two, order three and order four linear 

differential operators occurring in the suceptibility of the Ising model., (2010) 



Ising model, modular form,s and Calabi- Yau 44 

|htt p://www.uiiilini.fr/pagesperso/jacq ues-arthur.weil/Isingsolutions.htinl| 
[44] Y. Andre, Arithmetic Gevrey series and transcendence. A survey, Journal de Theorie des 

Nombres de Bordeaux, 15,(2003), 1-10. 
[45] Y. Andre, G-functions and geometry, Aspects of Mathematics E, Num. 013, ed. F. Vieweg in 

Braunschweig/Wiesbaden, (1989), ISSN: 0179-2156. 
[46] G. E. Andrews and B. C. Berndt, Ramanujan's Lost Notebook Part I (2005) and Part II (2009), 

Springer- Verlag, Berlin, New- York 
[47] N. D. Elkies, Explicit Modular Towers, Proceedings of the Thirty-fifth annual Allerton conference 

on communication, control and computing, ed. Basar, Tamer, and Alexander Vardy, 23-32. 

Univ. of Illinois at Urbana-Champaign 199, http://dash. harvard. edu/handle/1/2963974 
[48] W. P. Reinhardt, P. L. Walker, (2010), Weierstrass Elliptic and Modular Functions, NIST 

Handbook of Mathematical Functions, eds: F. W. J. Olver, D. M Lozier, R. F Boisvert and 

C.W Clark, Cambridge: Cambridge University Press 
[49] Y. Yang, On differential equations satisfied by modular forms, Mathematishe Zeitschrift 246, 

1-19 (2004). 
[50] Y. Yang and N. Yui, Differential equations satisfied by Modular Forms and K3 surfaces, Illinois 

Journal of Mathematics 51, Number 2, Summer (2007), 667-696 
[51] F. Beukers, Irrationality proofs using modular forms, Journees arithmetiques de Besangon, 

Asterisque No 147-148, (1987) 271-283, 345. 
[52] L. Long, Shioda-Inose structures of one-parameter families of K3 surfaces, J. Number Theory 

Journal of Number Theory 109, Issue 2, (2004), 299-318 
[53] F. Beukers and C. A. M. Peters, A family of K3 surfaces and C(3), J. Reine Angew. Math. 351, 

pp. 42-54, (1984) 
[54] J. Harnad amd J. McKay, Modular Solutions to Equations of Generalized Halphen 

type, R.Soc. Lond. Proc. Ser. A 456, Math. Phys. Eng. Sci, (2000), 261-294 and 

http://a rxiv.org/pdf/solv-int/9804006v3, 
[55] G. Halphen, Sur un systeme d'equations differentielles, C. R. Acad. Sci., Paris 92 (1881), 1101- 

1103. 
[56] R. Ivanov and M.P. Tuite, Rational Generalised Moonshine from abelian orbifoldings of the 

Moonshine Module, Nuclear Physics B 635, Issue 3, (2002) 435-472 
[57] M.P. Tuite, Monstuous Moonshine from Orbifolds, Comm. Math. Phys. 146, (1992) 227-309 
[58] G. E. Collins and M. J. Encarnacion, Efficient rational number reconstruction, J. Symbolic 

Comput. 20 (1995), no. 3, 287-297. 
[59] Y. Andre, Sur la conjecture des p-courbures de Grothendieck- Katz et un probleme de Dwork, 

Geometric aspects of Dwork theory. Editors A. Adolphson, F. Baldassarri, P. Berthelot, N. 

Katz and F. Loeser, Vol. I, II, 55-112, Walter de Gruyter, Berlin, New- York, (2004). 
[60] F. Beukers and G. Heckmann, Monodromy for the hypergeometric function „Fn—l, Invent. 

Math. 95 (2), pp. 325-354 (1989) 
[61] M. van der Put, M.F. Singer, (2003), Galois theory of linear differential equations, Grundlehren 

der Mathematischcn Wissenschaften, vol. 328, (Berlin: Springer) 
[62] G. Almkvist and W. Zudilin, Differential equations, mirror maps and zeta values, AMS IP 

Studies in Advanced Mathematics (2006), 38, 481-515 and arXiv:math/0402386 vr2 [math.NT] 
[63] B. Dwork, Differential Operators with Nilpotent p-Gurvature, American Journal of Mathematics, 

112, (1990) 749-786. 
[64] I. Bouw and M. Moeller Differential equations associated with nonarithmetic Fuchsian groups, 

J. London Math. Soc. (2) 81 (2010) 65-90, 

and arXiv/math/0710.5277vl [math.AG] 
[65] M. Dettweiler and S. Reiter On globally nilpotent differential equations, Journal of Differential 

Equations 248, Issue 11, (2010), 2736-2745 and arXiv/math/0605383 
[66] K.A. Driver and A.D. Love, Products of hypergeometric functions and the zeros of 4F3 

polynomials, Numerical Algorithms 26, 1-9, (2001). 
[67] J. Borwein and P. Borwein, Pi and the AGM: A Study of the Analytic Number Theory and 

Computational Complexity, Canadian Mathematical Society Series Monograph and Advanced 

Texts, John Wiley, New- York 1987. 
[68] Y. Andre and F. Baldassarri. Geometric theory of G-functions, Proc. Conf. in. Arith. Geom., 

Cortona, Symp. Math. XXXVII, 1-22, (1997), Cambridge: Cambridge University Press 
[69] J. Hadamard, Essai sur I'etude des fonctions donnees par leurs developpements de Taylor, 

Journal de Liouville, J. Math. Pures et Appliques. (4) , 8 (1892) 101-186 
[70] N. Zenine, S. Boukraa, S. Hassani and J.-M. Maillard, Ising model susceptibility: the Fuchsian 

differential equation for x '^'"■'i H^ factorizat ion properties, J. Phys. A: Math. Gen. 38 

(2005) 4149-4173 and |arXiv:cond-mat/0502155| 



Ising model, modular form,s and Calabi- Yau 45 

[71] G. Almkvist, C. van Enckcvort, D. van Stratcn and W. Zudilin, Tables of Calabi-Yau equations, 

[^Xiv:math0507430 http://arxiv.org/pdf/math/0507430 
[72] Y. Yang and W. ZudiUn, On Spi modularity of Picard- Fuchs differential equations for Calabi- 
Yau threefolds, with Yifan Yang, Contemporary Mathematics 517, 381-413, and arXiv: 

0803.3322v2 [math-NT] (2008) 
[73] A. J. Guttmann, Lattice Green functions and Calabi-Yau differential equations, J. Phys. A 42: 

Mathematical and Theoretical, (2009) 232001 (6pp), Fast Track Communication. 
[74] A. J. Guttmann, Lattice Green functions in all dimensions, J. Phys. A 43: Mathematical and 

Theoretical, (2010) 305205 (26pp). 
[75] V. Batyrev and M. Kreuzer, Constructing new Calabi-Yau 3-folds and their mirrors via conifold 

transitions, http://arxiv.org/pdf/0802.3376 
[76] Y-H. Chen, Y. Yang and N. Yui, Monodromy of Picard-Fuchs differential equations for Calabi- 
Yau threefolds, Journal fiir die reine und angewandte Mathematik (Crelles Journal), 616, 

167-203, (2008) and arXiv: 0605675v3 [math.AG] 
[77] S. Boukraa, S. Hassani, J-M. Maillard, B. M. McCoy, W. Orrick, N. Zenine, Holon- 

omy of the Ising model form factors, J. Phys. A 40, 75-111 (2007), and 

|http://arxiv.org/abs/math-ph/0609074 
[78] A. J. Guttmann, Spanning trees, Lattice Green functions and Calabi-Yau equations, 

http://insti.physics.sunysb.edu/scgp/conf/IntegrableModels2010/talks/Guttmann.pdf 
[79] G. Christol and B. Dwork, Differential modules of bounded spectral norms, Contemp. Math. 133 

(1992), 39-58. 
[80] G. Christol, Diagonales de fractions rationnelles, Seminaire de Theorie des Nombres, (1986/87), 

65-90. Progr. Math., 75, Birkhauser Boston, Boston MA (1988). 
[81] Y. Andre, G-fonctions et transcendance, J. Crelle, 476, (1996), pp. 95-125. 
[82] P. Griffiths, Periods of rational integrals on algebraic manifolds, I, II, Amer. J. Math. 90, (1968) 

pp. 568-626, 805-865, 

P. Griffiths, Periods of rational integrals on algebraic manifolds, III, Publ. Math. IHES;, 

(1970) 125-180. 
[83] Higher order Schwarzian ODEs. 

http://algo.inria.fr/bostan/HigherOrderSchwarzian.html' 

http://www.unilim.fr/pagesperso/jacqucs-arthur.wcil/HigherOrde rSchwarzian.htmll 

http://www. math. fsu.edu/^hoeij/maple. html/Higher Order Schwarzian. html' 
[84] C. Voisin, Symetrie Miroir, Panoramas et Syntheses, 2, (1996), Societe Mathematique de France. 
[85] L.Katzarkov, M.Kontsevich and T.Pantev, Hodge-theoretic aspects of mirror symmetry, Proc. 

Sympos. Pure Math., 78, Amer. Math. Soc, Providence, RI, 2008, pages 87-174. and 

arXiv:math.AG/0806.0107vl (2008) 
[86] M. Gross, D. Huybrechts and D. Joyce, Calabi-Yau Manifolds and Related Geometries, Lectures 

at a Summer School in Nordfjordeid, Norway 2001, (2003) Springer Universitext, Berlin, 

Heidelberg, New- York. 
[87] A. Strominger, S-T. Yau and E. Zaslow, Mirror Symmetry is T-duality, Nucl. Phys. B 479, 

243-259 and arXiv: hep-th/9606040 
[88] K. Hori and C. Vafa, Mirror symmetry, (2000) arXiv: hep-th/0002222 
[89] V. V. Batyrev and D. van Straten, Generalized Hypergeometric Functions and rational Curves 

on Calabi-Yau Complete Intersections in Toric Varieties, Comm. Math. Phys. 168 (1995), 

no. 3, 493533 and arXiv:alg-geom/9307010vl 
[90] N. Katz and T. Oda On the differentiation of the De Rham cohomology classes with respect to 

parameters, J; Math. Kyoto Univ 8, (1968) 199-213. 
[91] J. Harnad, Picard-Fuchs Equations, Hauptmoduls and Integrable Systems, Chapter 8 (Pgs. 137- 

151) of fntegrability: The Seiberg-Witten and Witham Equations (Eds. H.W. Braden and 

I.M. Krichever, Gordon and Breach, Amsterdam (2000)). 
[92] V. Golyshev, Classification problems and mirror duality, Surveys in geometry and number 

theory: reports on contemporary Russian mathematics, 88121, London Math. Soc. Lecture 

Note Ser. 338, Cambridge Univ. Press, Cambridge, (2007) and arXiv: math. AG/0510287 i/l 
[93] L. Schlesinger, Uber cine Klasse von Differentialsystemen beliebiger Ordnung mit festen 

kritischen Punkten, Journal fiir reine und angew. Math., 141, (1912), pp. 96-145. 
[94] D. Morrison, Picard-Fuchs equations and mirror maps for hypersurfaces, Essays on mirror 

manifolds, 241-264, Int. Press, Hong Kong, (1992) and arXiv:hep-th/9111025vl 
[95] K. Yamada, Pair Correlation Function in the Ising Square Lattice, Prog. Theor. Phys. 76, (1986) 

pp. 602-612. 
[96] J. R. Reyes Martinez, Correlation Function for the Z-invariant Ising models, Physics Letters 

A 227 (1997) pp.203-208, arXiv: hepth/9609135v2 (1997) 



Ising model, modular form,s and Calabi- Yau 46 

[97] J. R. Reyes Martinez, Multi-spin correlation functions for the Z -Invariant Ising model, Physica 

A 256 (1998) pp. 463-484 
[98] Y. Manin, Sixth Painleve equation, universal elliptic curve, and mirror 0/P2, AMS Transl. (2) 

186 (1998), 131-151. 
[99] J-M. Maillard and C-M. Viallet, Quasi- Integrability of the sixteen vertex model , Phys.Lett. B 

281, (1992), pp. 315-319 
[100] D. R. Morrisson, The Geometry Underlying Mirror Symmetry, New trends in algebraic 

geometry (Warwick, 1996), London Math. Soc. Lecture Note Ser., 264, 283-310 (1999) 

Cambridge Univ. Press, Cambridge, and arXiv:alg-geom/9608006i'2 
[101] N. M. Katz, Exponential Sums and Differential Equations, Annals of Mathematical Studies, 

Princeton Univ. Press. 124, (1990). 
[102] H. McKean and V. Moll, Elliptic Curves, Function Theory, Geometry, Arithmetic, Cambridge 

Univ. Press. First Published 1997 
[103] S. Hosono, A. Klemm and S. Theisen, Lectures on Mirror Symmetry, pp. 235-280 in Integrable 

Models and Strings, Proceedings, Helsinki, Finland 1993, Lecture Notes in Physics, 436 (1994) 
[104] M. van Hoeij, J-A. Weil, Solving second order linear differential equations with Klein's theorem. 

Proceedings of the 2005 International Symposium on Symbolic and Algebraic Computation 

(Beijing), 340-347, ACM, New York, 2005. 

[105] G. Moore, Attractors and arithmetic, arXiv:hep-th/9807056v2 (2003) 

[106] G. Moore, Strings and Arithmetic, Frontiers in Number Theory, Physics and Geometry II, 

Springer- Verlag Berlin Heidelberg (2007) cd. P. Cartier, B. Julia, P. Moussa, P Vanhove and 

arXiv:hep-th/0401049v4 (2007) 
[107] A. Strominger and E. Witten, New manifolds for superstring compactification, Comm. Math. 

Phys. 101 (1985), no. 3, 341-361. 
[108] E. Goursat, Memoire sur les fonctions hypergeometriques d'ordre superieur, Annales 

scientifiques de TENS, 2eme serie, tome 12, (1883), pp. 261-286 



